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ABSTRACT 


Lasers  have  been  incorporated  into  many  parts  of  society  including  heavy 
industrial  work,  military  defense,  light  displays,  scientific  research,  and  everyday 
household  appliances.  Free  Electron  Lasers  (FELs)  use  relativistic  electrons  traveling 
in  an  alternating  magnetic  field  to  produce  coherent  radiation.  One  particular 
application  for  a  FEL  is  to  beam  power  to  satellites.  The  United  States  and  Russia 
have  agreed  to  pursue  this  effort  together.  The  installation  in  the  United  States  is  the 
SpacE  Laser  ENErgy  (SELENE)  project  scheduled  for  installation  at  China  Lake, 
California. 

The  SELENE  laser  system  consists  of  an  oscillator  and  a  radiator.  The  electrons 
are  bunched  in  the  oscillator  in  preparation  for  entering  the  radiator.  The  use  of  a 
three-section  optical  klystron  has  been  proposed  for  the  SELENE  oscillator.  An  optical 
klystron  FEL  is  composed  of  multiple  undulator  sections,  each  separated  from  the  next 
by  a  dispersive  section.  This  dispersive  section  allows  for  larger  gain  from  a  shorter 
total  length  of  undulator.  The  effects  of  the  optical  klystrons  in  the  SELENE  oscillator 
using  the  proposed  SELENE  FEL  parameters  will  be  investigated  through  the  use  of 
simulations  of  two  and  three  section  undulators. 
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I.  INTRODUCTION 


The  United  States  and  the  Russian  governments  have  agreed  to  a  joint  venture 
to  investigate  the  possibilities  of  using  a  Free  Electron  Laser  (FEL)  to  provide  electrical 
power  to  satellites  from  a  ground  based  system  [1],  This  project  was  started  in 
Novosibirsk,  Russia.  The  initial  prototype  of  the  system  is  to  be  built  and  tested  at 
Novosibirsk  at  an  output  power  of  200  kW.  The  United  States  has  proposed  the 
SpacE  Laser  ENErgy  (SELENE)  system  to  be  built  at  China  Lake,  California.  This 
system  has  a  designed  output  of  1 0  MW. 

The  concepts  which  generated  the  need  for  the  SELENE  system  have  been 
around  since  the  early  days  of  the  space  program.  One  of  the  dominant  factors  is  the 
cost  per  kilogram  to  place  a  satellite  in  a  stable  orbit,  which  can  be  as  high  as 
$160,000  per  kilogram  for  a  geosynchronous  orbit  [1],  This  consideration  forces  the 
payload  to  be  mass  efficient.  Some  of  the  payload  must  be  used  for  the  power  and 
orbital  maintenance  systems.  The  power  systems  consist  of  large  solar  panels  used  to 
collect  sunlight  and  convert  it  to  electricity,  and  battery  backup  systems  used  to  store 
the  energy  for  use  when  the  satellite  is  in  a  blackout  period.  The  orbital  maintenance 
systems  generally  consist  of  chemically  powered  thrusters  used  to  maintain  a 
geosynchronous  orbit.  Due  to  the  limited  supply  of  chemical  fuel  for  the  thrusters,  the 
satellite- has  a  limited  time  to  maintain  orbit.  The  continuous  bombardment  of 
micrometeors  subject  the  solar  panels  and  connectors  to  physical  stresses  which  will 
render  them  inoperable.  The  continuous  deep  charge  and  discharge  cycle  of  the 
battery  system  will  also  limit  battery  lifetime.  Thus,  due  to  power  considerations, 
satellites  have  a  limited  lifetime,  considering  the  cost  to  place  them  into  orbit.  With  the 
large  number  of  satellites  and  satellite-supported  systems  in  the  military  today,  a 
means  of  reducing  the  cost  of  placing  satellites  into  orbit  and  extending  their  lifetimes 
would  be  of  great  benefit.  For  example,  it  would  allow  increased  functional  payload, 
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more  numerous,  smaller  satellites  for  greater  area  coverage,  or  significant  cost  savings 
which  could  be  used  for  technological  research. 

The  effect  SELENE  will  have  on  satellite  lifetimes  and  payload  is  tremendous, 
especially  if  a  system  of  similar  lasers  systems  is  installed  globally  to  provide  coverage 
for  all  satellites.  Continuous  coverage  for  all  geosynchronous  satellites  over  the 
United  States  may  be  achieved  with  three  stations.  With  the  use  of  a  beacon  satellite 
to  direct  the  beam  to  individual  satellites,  all  Low  Earth  Orbit  satellites  in  the  Western 
Hemisphere  north  of  the  equator  can  receive  power  from  SELENE  located  at  China 
Lake,  California.  [2] 

The  current  solar  cells  used  in  satellites  are  inefficient.  The  use  of  silicon 
photocells  is  common  in  most  satellites.  These  cells  receive  broad-band  solar  radiation 
and  convert  10-15%  of  the  incident  radiation  to  electricity  for  use  by  the  satellite  [2]. 
Using  the  SELENE  project  to  tune  the  light  to  a  wavelength  near  the  peak  sensitivity  of 
the  current  silicon  cells  can  improve  the  overall  conversion  efficiency  by  a  factor  of 
three.  The  capability  of  focusing  the  laser  on  the  solar  cells  will  increase  the  incident 
intensity  on  the  solar  cells  by  a  factor  of  4.6  allowing  for  smaller  arrays  of  photocells  to 
achieve  the  same  power  production.  Using  the  current  silicon  based  solar  cells,  this 
represents  a  gain  of  almost  14  in  power  output  per  unit  area  of  solar  panel.  With 
SELENE,  the  solar  cells  can  be  designed  for  an  optimum  wavelength  for  electrical 
conversion  improving  the  conversion  efficiency  by  a  factor  of  four.  These  smaller 
arrays  will  reduce  the  weight,  as  well  as  the  cross-sectional  area  of  the  satellite.  This 
will  reduce  the  number  of  interactions  the  satellite  has  with  micrometeors  which  will 
enhance  the  reliability  of  the  satellite’s  power  systems. 

The  chemical  propulsion  systems  required  for  station  keeping  may  also  be 
replaced  with  a  system  of  electrically  powered,  ion  thruster  propelled,  satellite  tugs. 
These  tugs  can  be  used  to  precisely  place  the  satellite  in  the  exact  orbit  required  for 
them  to  perform  their  mission.  This  will  remove  the  problem  of  fuel  supply  exhaustion 
as  well  as  removing  the  large  mass  required  by  the  chemical  boosters  and  storage 
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systems.  With  these  changes  in  standard  satellite  design,  the  mass  of  the  support 
systems  on  the  satellite  and  the  cost  per  kilogram  to  orbit  will  be  reduced.  The  return 
on  satellite  lifetime  is  large,  considering  the  expense  of  placing  a  single  satellite  in 
orbit  multiplied  by  the  enormous  number  of  satellites  available.  This  innovative  system 
will  have  a  cost  return  of  almost  ten  to  one.  [1]. 

The  SELENE  project  design  is  unique  in  its  three-section  undulator.  The 
purpose  of  the  dual  dispersive  sections  in  the  undulator  is  to  provide  additional 
flexibility  in  the  design  and  operation  of  the  FEL,  as  well  as  to  reduce  the  steady-state 
optical  field  strength  in  the  oscillator.  The  separate  oscillator  and  radiator  sections  are 
necessary  to  provide  bunching  of  the  electrons  in  a  weak  optical  field  prior  to  their 
entering  the  radiator.  This  allows  for  large  power  output  from  the  radiator,  which  does 
not  require  optical  devices,  while  keeping  the  power  density  low  in  the  oscillator  to 
minimize  damage  to  the  optical  devices  in  the  system.  The  beam  from  the  radiator 
then  travels  through  an  underground  vacuum  tube  where  it  expands  sufficiently  for  the 
mirrors  to  project  it  to  the  satellites  with  minimal  mirror  damage. 


3 


II.  FREE  ELECTRON  LASER  THEORY 


A.  INTRODUCTION 

The  basic  design  of  an  FEL  oscillator,  shown  in  Fig.  1 ,  consists  of  a  relativistic 
electron  beam,  a  wiggler,  and  an  oscillator  cavity.  The  wiggler,  or  undulator,  consists 
of  a  series  of  periodic  magnetic  fields  either  in  a  linear  configuration  or  in  a  helical 
pattern.  The  helical  undulator  is  more  efficient  at  transferring  energy  from  the 
electrons  to  the  optical  field,  but  the  linear  configuration  is  much  easier  and  less 
expensive  to  construct.  Most  operating  FEL’s  today  use  a  linear  undulator.  The 
oscillator  consists  of  two  mirrors,  one  of  which  is  completely  reflective  and  the  other  is 
partially  reflective  to  allow  optical  power  to  be  extracted. 

Due  to  the  transverse  acceleration  provided  by  the  magnetic  fields  in  the 
undulator,  the  electrons  in  the  FEL  emit  radiation  in  a  cone  directed  along  the 
electron’s  velocity  vector.  The  size  of  the  cone  is  related  to  the  velocity  of  the 
electrons.  Some  of  this  radiation  is  trapped  by  the  oscillator  and  provides  the  optical 
field  for  interaction  with  following  electrons.  The  interaction  between  the  electrons,  the 
magnetic  field,  and  the  optical  field  cause  bunching  in  the  electron  beam  on  the  scale 
of  the  optical  wavelength  of  the  FEL.  This  microscopic  bunching  produces  coherent 
radiation  and  gain  in  the  optical  field.  This  is  the  basic  principle  behind  the  FEL. 

With  the  large  number  of  high  energy  electrons  present  in  the  electron  beam  and 
the  relatively  large  acceleration  produced  by  the  undulator  fields,  the  possibility  of  high 
output  power  is  easily  recognized.  The  frequency  of  the  radiated  light  is  proportional 
to  the  energy  of  the  accelerating  electrons,  thus  providing  for  coherent  light  which  can 
be  tuned  over  a  wide  range  of  frequencies.  Some  of  the  operating  parameters  of  an 
FEL  are  the  number,  size  and  magnetic  field  strength  of  the  undulator  periods,  the 
electron  beam  energy,  current,  and  pulse  shape,  the  quality  of  the  electron  beam,  the 
transmission  coefficient  of  the  oscillator  output  mirror,  the  presence  of  special 
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undulator  designs  such  as  dispersive  sections,  length  of  the  oscillator,  and  the  radius 
of  the  mirrors.  The  effects  these  parameters  have  on  FEL  design  and  output 
characteristics  will  be  studied,  with  emphasis  on  the  effects  of  dispersive  sections  on 
high  gain  FELs. 


ELECTRON  ELECTRON 

ACCELERATOR  BENDING  MAGNETS  BEAM  DUMP 


OUTPUT 


Figure  1 .  Elements  of  a  typical  FEL  oscillator. 


B.  RADIATION  FROM  CHARGED  PARTICLES 

For  a  single  electron,  the  Lienard-Wiechert  potentials  in  SI  units  are  [3], 


gjfp  %r) 

471 


9^r)  . 


(2.2) 


where  §LW  is  the  electric  field  potential,  ftLW  is  the  vector  potential,  e  is  the  charge  on 
the  electron,  ~S{f,xr)  is  the  vector  from  the  source  point  to  the  field  point  T,  xr  is  the 
retarded  time  at  which  the  particle  radiated  in  order  for  the  radiation  to  reach  t  at  time 
t,  V(xr)  is  the  velocity  of  the  electron  at  time  xr,  e0  is  the  permittivity  of  free  space,  (Xq 
is  the  permeability  of  free  space,  and 
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(2.3) 


=  [1-^(7^)  ^(t,)]-1 
c 

is  the  relativistic  correction  factor  where  c  is  the  speed  of  light  and  7t  is  the  unit  vector 
in  the  direction  of  3".  The  electric  (£*)  and  magnetic  (£?  )  fields  can  be  derived  from 
Maxwell’s  equations, 

=  Vx/S LW(t,t)  (2,4) 

and 

£ (T,f)  =  -V<| )LW  ~^-XLW  .  (2.5) 

Each  field  can  be  broken  down  into  two  components  as  shown  below  in  Equations  2.6 
and  2.7,  one  dependent  on  the  velocity  of  the  particle  and  the  other  proportional  to  the 
acceleration  as  denoted  by  the  respective  subscripts  v  and  a.  These  field  vectors  are 
[3] 

£,W)  -  ~  4  c2 

471  sz 

4ji  S  C 

i  LL  J 

and 


= 

c 

^(7,0  =  Bv  +  £a  '  .  (2.7) 

=  ^~7tx£a 
c 

The  electric  field  proportional  to  the  acceleration,  Ea,  has  a  directional 
component  f  (0,v),  which  can  be  defined  as 
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(2.8) 


f(Q,v)s 


rt-- 


1  +-C-2— cos0 


1/2 


1 - COS0 

c 


where  0  is  the  angle  between  V  and  7t.  Looking  in  the  direction  of  V  where  0  =  0, 
f  (0,v)  =  (l-v/c)'2  which  becomes  very  large  as  v  approaches  c.  Thus,  in  the 
extreme  relativistic  limit,  the  radiation  will  normally  be  peaked  in  the  direction  of  V. 
This  is  not  always  the  case  as  there  are  additional  factors  not  noted  above. 
Specifically,  in  the  case  that  the  acceleration  is  in  the  same  direction  as  the  velocity, 
(ft  -  V/c)  x  ^  =  0  for  a  field  point  in  the  direction  of  V,  and  there  is  no  radiation  emitted 
in  the  direction  of  V.  For  the  FEL,  the  undulator’s  magnetic  field  produces  a  force  on 
the  electron  which  is  perpendicular  to  the  motion  of  the  electron  and  /(0,t?)  is 
sufficient  to  describe  the  radiation  pattern  produced  by  a  single  electron. 

Table  1  shows  some  typical  ranges  for  the  angular  spread  of  the  radiation 
pattern,  in  degrees,  where  0'(x)  is  the  angle  where  f{d',v)  =  xf{ 0,v),  Te  is  the 
relativistic  kinetic  energy  of  the  electron  in  MeV,  and  x  is  the  ratio  of  the  directional 
factor  at  the  angle  0'  to  the  maximum  value  of  the  directional  factor  at  0  =  0.  For 
example,  for  an  electron  with  y=100,  the  directional  factor  /(0',v)  falls  to  90%  of 
f  (0,v)  at  an  angle  of  0'  =  2.5°  from  V  as  shown  in  bold  type  in  Table  1.  The  Lorentz 
factor,  y=(l-^2)  1/2  where  =  v7c,  will  be  used  extensively  in  the  development  of 
FEL  theory  since  the  electron  beam  is  relativistic. 
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V/C 

y 

Te 

(MeV) 

O'  (0.5) 

(degrees) 

Q'  (0.1) 

(degrees) 

6'  (0.01) 

(degrees) 

0.90000 

2.29 

.66 

29.1 

51.1 

98.2 

0.99000 

7.09 

3.1 

13.1 

20.1 

34.1 

0.99900 

22.37 

10.9 

5.6 

8.2 

13.5 

0.99990 

70.71 

35.6 

2.3 

3.3 

5.4 

0.99995 

100.0 

50.6 

1.79 

2.5 

4.1 

0.99999 

223.6 

114 

.96 

1.3 

2.2 

Table  1 .  Angular  spread  of  the  radiation  from  a  relativistic  electron.  After  Ref.  [3], 


C.  ELECTRON  INTERACTION  WITH  MAGNETIC  FIELDS 

The  path  taken  by  the  electrons  while  in  the  undulator  determines  the  FEL 
characteristics.  The  electrons  are  injected  into  the  undulator  at  relativistic  velocities 
and  encounter  the  magnetic  fields  of  the  undulator,  producing  spontaneous  radiation. 

The  magnetic  field  of  the  undulator  is  usually  produced  by  permanent  magnets  or 
current  carrying  coils.  The  electrons  are  directed  along  the  axis  of  the  undulator, 
which  will  be  represented  as  the  z-axis.  For  linear  undulators,  the  magnets  will  be 
aligned  with  the  x-axis.  When  the  electrons  are  injected  into  the  undulator,  they  see 
the  magnetic  field  as  being  Lorentz-contracted.  In  a  helical  undulator  with  circular 
polarization,  the  magnetic  field  in  cylindrical  coordinates  ( r,Q,z )  is  [4] 
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BrH=2B0 


sin(0-/coz)  , 


IqM 


MV)' 

kQr 


B 


H 

e 


2B0 

V 


li{k0r)cos(Q-k0z),  and 


(2.9) 


SZH  =  — 2B0/1(/c0r)  cos(0-/coz) , 

where  S0  is  the  peak  magnetic  field  strength  of  the  undulator,  kQ  =  2ttA0  is  the  wave 
number  of  the  undulator,  and  is  the  length  of  one  undulator  period.  /0  and  /1  are 
hyperbolic  Bessel  functions  [5].  Off  the  axis,  it  is  the  transverse  components  of  Eq. 
2.9  which  act  to  modify  the  electrons  trajectory  along  the  undulator  length.  Near  the 
axis,  Eq.  2.9  may  be  expanded  for  k0r<£  1  to  find  the  increase  in  the  average 
transverse  magnetic  field  as  r  increases, 

B h  ~  Bq(1 +/C()f2/4  +  •  •  •  )  .  (2.10) 

This  higher  magnetic  field  density  away  from  the  axis  results  in  the  electrons  being 
tunneled  back  towards  the  axis  where  the  ideal  on-axis  magnetic  field  in  rectangular 
coordinates  is 

=  [Bx,By,Bz]  =  S0[cos(/c0z),sin(/r0z),0]  .  (2.1 1) 

The  magnitude  of  the  transverse  magnetic  field  in  an  helical  undulator,  is  constant 
along  the  z  axis,  =  B0.  This  allows  for  simpler  analytical  analysis  of  FEL 
interactions. 

For  a  linear  undulator,  the  on-axis  magnetic  field  in  rectangular  coordinates  [4]  is 

EL  =  jB0[0,sin(/r0z)cosh(/c0y),cos(/r0z)sinh(k0y)]  (2.12) 

and  the  increase  in  the  average  magnetic  field  strength  as  the  electron  moves  off  axis 
is 


B, 


BL  -  -^(1  +k02y2/2+  •  -  • ). 


(2.13) 
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As  can  be  seen  by  Eq.  2.13,  the  linear  undulator  provides  for  focusing  in  the  y 
direction  only.  Radial  focusing  may  be  achieved  using  quadrupole  lenses  or  by  using 
magnetic  poles  with  parabolic  faces  instead  of  flat  faces  [6].  The  ideal  magnetic  field 
strength  on  axis  for  a  linear  undulator  is 

=  B0[0,sin(/c0z),0],  (2.14) 

and  the  average  magnetic  field  strength  seen  by  the  electrons  is  El  =  B q/^2.  This 
results  in  smaller  acceleration  in  the  transverse  direction  as  compared  with  the  helical 
undulator  and  thus  reduced  emissions,  coupling  and  gain.  The  helical  undulator 
naturally  provides  radial  focusing  without  the  added  engineering  required  by  the  linear 
undulator.  From  the  evidence  presented,  it  would  seem  that  the  helical  undulator 
design  is  superior  in  performance  as  compared  to  the  linear  undulator.  However,  the 
helical  undulator  uses  current  carrying  coils  for  creating  the  magnetic  field  in  the 
undulator  whereas  the  linear  undulator  uses  permanent  magnets.  Thus  the  cost  of 
constructing  and  operating  helical  undufators  results  in  the  use  of  linear  undulators  in 
most  modern  FELs. 

For  an  electron  perfectly  injected  into  an  ideal  linear  undulator,  using  the  Lorentz 
force  equation  in  CGS  units, 


q(£+txB), 


(2.15) 


the  magnetic  field  of  (2.14),  and  assuming  y  «  constant,  results  in 


P*  =  Pj.  =  Pz  B0sin(/r0z) 


Py  =  0 


(2.16) 


P 


Z 


q 

mcy 


P*  B0s\n(k0z) 


for  equations  of  motion  in  the  absence  of  an  optical  field,  where  m  is  the  mass  of  the 
electron  and  q  =  -  e  is  the  charge  on  the  electron.  Using  z  =  &zct  and  the 


11 


approximation  (3Z  ~  constant  ~  1  together,  it  can  be  shown  that 


x  = 


-—sin (k0ct),  and 
ykQ 


X  =  PjC  = 


Kc 

y 


cos(k0ct)  . 


(2.17) 


This  gives  a  peak  displacement  in  the  x  direction  of  K/k0 y,  where  K  =  eB0Xol2nmcz  is 
the  dimensionless  FEL  parameter.  For  a  typical  FEL,  y  =  100,  e=  4.8x1 0-10 
statcoulombs,  m=9.1x10-28  gm,  pz»l,  c  =  3x1010cm/s,  X0  =  5  cm, 
k0  =  2x1X0  =  1.25cm"1,  and  B0  =  4  kG  ,  which  results  in  K=1.87  and  a  peak 
displacement  in  the  x  direction  of  0.15  mm. 


In  order  for  the  proper  coupling  between  the  electron  beam  and  the  optical  field 
to  occur,  the  electron  beam  should  stay  within  the  optical  mode  while  in  the  undulator. 
From  a  basic  study  of  optics,  the  distance  over  which  an  optical  mode  area  doubles 
due  to  natural  diffraction,  called  the  Rayleigh  length,  is  z0  =  nw§ A,  where  w0  is  the 
radius  of  the  optical  field  at  the  waist.  If  the  Rayleigh  length  is  not  matched  to  the 
undulator,  the  optimal  coupling  cannot  occur.  For  example,  if  the  Rayleigh  length  is 
much  smaller  than  the  undulator  length,  then  the  optical  mode  area  will  become  much 
larger  than  the  electron  beam.  This  will  result  in  reduced  coupling  for  two  reasons. 
First,  the  optical  field  strength  will  be  reduced  in  amplitude  since  it  will  be  spread  over 
a  larger  optical  mode  area  and  second,  the  electrons  will  be  interacting  with  a  smaller 
fraction  of  the  total  mode  area.  Thus,  the  capability  of  the  undulator  and  optical  field 
to  bunch  the  electrons  and  produce  coherent  radiation  will  be  reduced.  For  a  Rayleigh 
length  equal  to  the  length  of  the  undulator,  the  optical  mode  area  will  be  *  XL,  which, 
for  a  typical  FEL  with  X  -  1 0  p.m  and  L  =  5  m  results  in  a  mode  radius  of  4.0  mm. 
For  a  typical  FEL,  the  wiggle  motion  of  the  electrons  in  the  undulator  remains  inside 
the  optical  mode.  The  relative  size  of  the  electron  beam  with  respect  to  the  optical 
field  is  quantified  by  the  filling  factor  F  -  nrblnwQ  where  rb  is  the  radius  of  the 
electron  beam. 
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D.  ELECTRON  INTERACTION  WITH  OPTICAL  FIELDS 

The  interaction  of  the  electrons  with  the  light  field  determines  the  FEL’s  ability  to 
produce  coherent  light.  As  the  electrons  travel  through  the  undulator,  they  oscillate  in 
the  transverse  direction  once  per  undulator  period.  The  electrons  must  move  slower 
than  the  light  in  the  undulator.  When  the  change  in  relative  position  of  the  electrons 
and  optical  field  is  equal  to  one  wavelength  of  light  at  the  end  of  one  undulator  period, 
the  electrons  are  at  resonance  [8].  Therefore  X=c(1-p2)Af  where  X  is  the 
wavelength  of  the  light  emitted  at  resonance,  and  At  =  X0lc&z.  Solving  for  p2  using 
the  Lorentz  factor,  and  pf  =  K2/^,  results  in  p2  *  1  -(1  +  K2)l2f.  The  resonance 
condition  follows  directly, 


1  1  +  K2  .  . 

*-**-&-  (218) 

This  shows  that  the  output  frequency  of  the  FEL  can  be  tuned  by  changing  the 
undulator  parameters  X0  or  K,  or  more  readily,  by  changing  the  electron  beam  energy, 
ymc2. 

The  motion  of  the  electrons  in  the  undulator  may  be  described  by  a  simple 
pendulum  equation  as  shown  below.  The  plane  electromagnetic  wave  propagating 
down  a  helical  undulator  may  be  described  by  [4] 

tr  =  E(z,t)  [cos(\|/),~ sin(\|/)  ,0] 

&r  =  =  E{z,t)  [sin(\|/),cos(\y),0]  .  (2‘19) 

where  \j/  =  kz-<ot+§,  k  =  2it/X,  co  =  kc  is  the  angular  frequency  of  the  light,  X  is  the 
optical  wavelength,  <j>  is  the  phase  of  the  light,  and  E  is  the  magnitude  of  the  optical 
wave.  The  behavior  of  the  relativistic  electrons  in  the  presence  of  electromagnetic 
fields  can  be  described  using  the  Lorentz  force  equation,  (Eq.  2.15)  and  the  time  rate 
of  change  of  the  energy  of  the  electron, 


^®2fL  =  -eci?-£, 


(2.20) 
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where  the  energy  of  the  electron  is  y me2.  In  total,  the  electron  interacts  with  the  fields 
^  ~  and  &  =  &0  +~§r.  Solving  the  Lorentz  force  equation  for  assuming  perfect 
injection  into  a  helical  undulator  and  p2  ~  1  results  in 

Pi  =  -  —[COS  (/roz),sin(V).0]  •  (2.21) 

Thus,  inserting  Eq.  2.21  and  tr  from  Eq.  2.19  into  Eq.  2.20  gives 

Y=f^C0S(C+<t’)'  (2.22) 

where  the  relative  phase  between  the  electron  and  the  optical  field  is  given  by 
C  =  (k  +  k0)z -at.  For  -n/2  <£+<}><  n/2,  y>0  and  energy  is  transferred  from  the 
optical  field  to  the  electrons,  and  for  n/2  <  £+<j>  <  3n/2,  y<  0  and  energy  is  transferred 

from  the  electron  beam  to  the  optical  field.  Now  the  effect  of  bunching  becomes 

apparent.  If  the  electrons  are  bunched  near  £+4 >=7t,  then  there  will  be  a  large 
fraction  of  electrons  radiating  energy  to  the  optical  field  in  the  same  phase  and  the 
optical  field  will  grow.  If  the  same  group  of  electrons  are  bunched  near  £+<j>  =  0,  then 
the  electrons  will  be  absorbing  energy  from  the  optical  field  and  the  optical  field 
strength  will  decrease.  If  the  electrons  are  evenly  distributed  in  phase  with  respect  to 
the  light,  then  as  many  electrons  will  lose  energy  as  will  take  energy  from  the  optical 
field  and  the  net  change  in  optical  field  strength  will  be  zero. 

Using  the  definition  of  y,  Eq.  2.17,  differentiating  y2  with  respect  to  t,  and 
equating-  with  Eq.  2.22,  the  fractional  change  per  unit  time  in  the  energy  of  the  electron 
is  found  to  be 


X  = 

Y  1  +K 2 


eEK  „  tY  ,, 

-5 - COS(C  +  <|>)  • 

- fmc 


(2.23) 


Taking  the  second  derivative  of  £  with  respect  to  time,  solving  for  p2,  and  substituting 
the  result  into  Eq.  2.23  gives 


(1  +Kz){k+kQ)c 


t  = 


-f^cos(C+<i>)  - 


(2.24) 
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Using  the  resonance  condition  (Eq.  2.18),  ~  1  and  k  »  k0  for  y  »  1,  the  pendulum 

equation  is  found  to  be 


C  * 


2k0eEK 

■fm 


cos(t;+<j>)  . 


(2.25) 


In  order  to  simplify  calculations  and  to  get  a  scale  that  is  valuable  to  compare 
different  FEL  designs,  the  use  of  dimensionless  parameters  is  introduced.  A 
dimensionless  time,  x,  is  defined  by 


(2.26) 


where  L  =  N\q  is  the  length  of  the  undulator,  N  is  the  number  of  undulator  periods 
and 


w 2  r2  H2  r2  00 

-=—(  )  =  — — —  (  )  -  —(  ) 

dt2{  f  L2d^-]  L2  ” 


(2.27) 


This  "normalizes"  the  FEL  so  that  x  goes  from  0  to  1  as  the  electrons  travel  down  the 
undulator.  The  dimensionless  pendulum  equation  for  a  helical  undulator  then 
becomes 


V  0  AnNeKLE  ,r  ,  ,  ,Y 
C  =  v  =  cos(C+<|))  =  lalcos(C+<j))  (2.28) 

where  la  I  is  the  magnitude  of  the  dimensionless  optical  field  strength  in  the  oscillator 

o 

and  v  =  C  is  the  relative  phase  velocity  between  the  optical  field  and  the  electron. 

For  a  linear  undulator,  the  electromagnetic  radiation  in  the  undulator  is  of  the 
form  [4], 


=  E(z,f)[cos(/(z-cof  +<|>(z,f)),0,0] 
4-C*\0  =  E (z,f)[0,cos(/cz  -cof  +<|>(z,f)),0] . 


(2.29) 


This  radiation  field,  combined  with  the  magnetic  field  of  the  linear  undulator  from  Eq. 
2.12,  will  result  in  a  similar  dimensionless  pendulum  equation. 
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where  J0  and  J1  are  solutions  to  Bessel’s  differential  equations,  £  =  /C2/  2(1  +  K2),  and 
K  =  eB0X0l'l22nmcz  is  the  rms  undulator  parameter. 


E.  OPTICAL  WAVE  DEVELOPMENT 

To  see  how  the  optical  fields  develop  and  interact  with  the  electron  beam,  the 
derivation  starts  with  Maxwell’s  equations  in  CGS  units.  The  relationship  between  the 
vector  potential,  and  the  current  density  perpendicular  to  the  direction  of 

propagation  of  the  optical  field,  Jjf.t)  as  functions  of  position,  t  and  time,  t,  can  be 
used  to  derive 

J-?  ■  (2.3i) 

The  vector  potential  for  a  circularly-polarized  plane  wave  propagating  in  the  1  direction 
is  given  by  t  =  (1//c)  E(z,t )[sin(\y)  ,cos(\j/)  ,0],  where  v| /  =  {kz-a>t +<Kz,f))  and  ©  =  kc 
is  the  angular  frequency  of  the  optical  field.  The  operator,  ^2,  may  be  separated  into 
transverse  and  longitudinal  components,  =  ^z+tfz2.  For  initial  development  only 
perfect  injection  and  plane  wave  propagation  will  be  used,  for  which  ^2  =  Vz2  applies. 
It  is  assumed  that  the  optical  field  envelope  evolves  slowly  over  time  with  respect  to 
the  optical  period  and  varies  slowly  over  space  with  respect  to  the  optical  wavelength. 
This  is  called  the  slowly-varying  amplitude  and  phase  approximation  and  is 
summarized  by 

BE  r-  8<j) 

IF*”*  '  IF *“*  ' 


(2.32) 


To  solve  Eq.  2.31,  each  component  of  the  left-hand  side  is  evaluated.  The 
results  and  the  right-hand  side  are  projected  onto  perpendicular  vectors  which  will 
allow  separation  of  the  solution  to  Eq.  2.31  into  two  separate,  but  dependent 
equations.  The  first  step  is  to  take  the  second  partial  derivative  of  X  with  respect  to  z 
which  results  in 


0“  fc+H 


+  I  If +f  0  • 


(2.33) 


Next,  the  second  partial  derivative  of  X  with  respect  to  t  is  divided  by  c2  to  get 


1dX_  1  a2E  E 


2  dt  c2k  dt2  kc2  dt 


2  BE  jty  _  E  92(j> 

c2k  dt  dt  c2k  dt2 


CD  [sin(\|/) ,  cos(vy) ,  0] 


[cos(\y) ,  -sin  (y) ,  0] 


(2.34) 


Adding  Eq.  2.33  and  Eq.  2.34  together  and  applying  Eq.  2.32,  results  in  all  of  the 
second-order  terms  dropping  out  and  the  resulting  equations  reduce  to 


~2E  H  -  !  If  rinM '  °osM .  0] 


+  2  dz  +  "c  IF  *-C0S^  ’  ~S'n  ^ 1  • 


(2.35) 


If  two  unit  vectors  are  defined  by  e,  =  [cos(\y)  ,-sin(\ v) ,  0]  and  £2  =  [sin(y) ,  cos(\)/) ,  0], 
Eq.  2.35  is  combined  with  the  right-hand  side  of  Eq.  2.31 ,  and  the  dot  products  of  e, 
and  £2  are  found  for  both  sides  of  the  resulting  equation,  then  the  result  is 
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(2.36) 


_  jl_ 

dz  C  dt 


hj,- 


i<£2 


The  current  density  is  composed  of  the  sum  of  the  individual  electron  currents. 
Combining  this  sum  with  from  Eq.  2.21 ,  the  right  hand  sides  of  Eq.  2.36  become 


^.•£i  =  ^—'Lcos(\y+k0zi)b{3)(f-Ti) 

y  i 

^  =  -^^£sin(\|/+/c0Z;)5(3)(T-7}) 

y  i 


(2.37) 


where  7/  is  the  position  of  the  Ith  electron.  To  get  a  useful  approximation  for  the 
current  density,  an  average  over  all  of  the  electrons  weighted  by  the  number  density  of 
the  electrons  in  the  electron  beam  p  is  used.  If  it  is  assumed  that  a  long  electron 
pulse  relative  to  the  wavelength  of  the  light  is  used  and  that  there  are  no  optical 
modes  along  the  length  of  the  undulator,  then  the  approximations  d£/3z-» 0  and 
dtydz ->0  can  be  made.  These  approximations  and  the  averaging  of  7lt  along  with  the 
definitions  of  vy  and  t,  are  used  to  find 


ox  y 

E&-  =  2l»^<sin(C+«> 

ch  y 

which  together,  are  the  real  and  imaginary  parts  of 

JL  Eg1*  =  2neKL  p  <c-/T> 

0T  Y 


(2.38) 


(2.39) 


If  a,  the  dimensionless  complex  optical  field  strength,  and  j,  the  dimensionless  current 
density,  are  defined  by 


„  AnNeKLEe * 
a  =  - ^ - t - 

-fmc2 

(2.40) 

■  8N(enKL)2p 
J  fmc2  ’ 

then  Eq.  2.38  can  be  rewritten  as 

a=-j<e-^>,  (2.41) 

which  is  known  as  the  dimensionless  FEL  wave  equation  and  describes  the  evolution 
of  the  optical  wave  as  it  propagates  along  the  undulator  in  phase  space.  The  results 
are  similar  for  a  linear  undulator  [4],  where  the  FEL  parameter  K  is  replaced  by 
K(J0^)-J^))  and  E  is  replaced  by  E  =  EH 2  in  the  terms  a  and  j  in  the  preceding 
derivation.  For  typical  FELs,  the  efficiency  is  low,  on  the  order  of  1 0%. 

F.  ELECTRON  BEAM  QUALITY 

For  electrons  not  injected  perfectly,  the  focusing  of  the  electrons  back  towards 
the  axis  produces  betatron  oscillations  [4],  The  harmonic  betatron  oscillations  have  a 
constant  of  motion  [4],  =  [(ogy2(i)+y2(i)]/2  =  [cop2yo  +  y2(0)]/2  where  y0  is  the  initial 

electron  transverse  position,  y(0)  =  L0y,  and  0y  is  the  initial  injection  angle.  An 
electron  injected  off-axis  at  y0  or  at  an  angle  of  0y  will  experience  an  average  phase 
velocity  change  of 

Av  a  =  -JZtL(K2k$y$  +t202)  (2.42) 

For  lAvpl  >  re,  the  gain  bandwidth,  the  FEL’s  ability  to  bunch  the  electrons  and 
produce  gain  will  become  inhibited. 

The  measure  of  the  ability  of  an  electron  accelerator  to  inject  electrons  into  an 
FEL  with  the  proper  initial  position  and  velocity  is  called  emittance,  sy  =  yo0o  for  a 
linear  undulator  and  e  =  rbQ0  for  a  helical  undulator  where  y0  and  rb  are  the  rms  initial 
positions  and  0O  is  the  rms  initial  angular  deflection  of  the  electrons  away  from  the  z 
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axis.  For  any  accelerator,  either  y0  or  0O  may  be  reduced  by  external  focusing,  but 
the  product  e  is  fixed.  In  order  to  minimize  Av  in  Eq.  2.42,  the  angular  and  position 
spreads  are  matched  so  that  Kk0y0  =  y6y.  For  a  beam  with  matched  emittances,  Eq. 
2.42,  and  the  resonance  condition,  Eq.  2.18,  can  be  used  to  show  that 
Avp  ~  2nNKkey  / y.  Since  it  is  required  that  lAvpl  <  n  in  order  to  maintain  proper  FEL 
coupling,  then  e™ax  ~  y/NKk  =  yk/2nNK.  Typical  FEL’s  require  a  maximum  emittance 
of  1  mm-mradian.  With  this  limitation  on  electron  beam  quality,  the  effectiveness  of 
the  FEL  at  producing  coherent  light  is  dependent  on  several  design  parameters. 
These  parameters  are  more  restrictive  when  the  optical  wavelength,  k,  is  short  or 
when  the  number  of  undulator  periods,  N,  is  large. 

Inadequate  beam  quality  can  degrade  FEL  performance  by  inhibiting  electron 
bunching.  For  the  studies  on  SELENE,  the  electron  beam  is  considered  to  be  injected 
into  the  oscillator  with  initial  electron  energies  distributed  about  v0  in  a  Gaussian 
distribution  with  standard  deviation  of  oG  =  4n/VAy/y.  The  injection  angle  is  assumed 
to  be  an  exponential  distribution  centered  at  zero  with  a  standard  deviation  of 
cre  =  4Tt/Vy202/  (1  +  K2). 

G.  PHASE-SPACE  AND  THE  LOW  CURRENT,  LOW  GAIN  FEL 

Since  the  pendulum  equation  is  periodic  within  an  optical  wavelength,  the 
evolution  of  the  electron  with  respect  to  its  phase  and  phase  velocity  relative  to  the 
optical  field  over  a  distance  of  one  wavelength  of  light  can  be  studied  in  order  to 
determine  the  overall  characteristics  of  the  FEL.  Figure  2  shows  the  evolution  of  the 
electrons  as  they  traverse  the  undulator.  The  electrons  are  initially  injected  into  the 
undulator  on  resonance,  v0  =  0.  From  Eq.  2.23,  electrons,  whose  phase  with  respect 
to  the  optical  field  is  -n/2<l<nl2,  will  accept  energy  from  the  optical  field.  These 
electrons  are  depicted  on  the  left  side  of  the  v  versus  £,  or  phase-space,  plot.  The 
electrons  are  shaded  so  that  at  the  beginning  of  the  undulator  they  are  a  light  grey, 
and  become  darker  as  they  travel  down  the  undulator.  Thus,  as  the  electrons  on  the 


20 


left-hand  side  of  ine  phase-space  plot  travel  down  the  undulator,  they  gain  energy. 
This  increases  their  phase  velocity  and  they  begin  to  advance  in  phase  with  respect  to 
the  optical  field.  For  the  electrons  on  the  right-hand  side  of  the  phase  space  plot, 
tc/2< £<3tt/2.  They  lose  energy  to  the  optical  field  as  they  propagate  down  the 
undulator.  Due  to  this  decrease  in  energy,  their  phase  velocity  is  reduced  and  they 
retard  in  phase  with  respect  to  the  optical  field.  These  effects  combine  to  produce  the 
clockwise  rotational  effect  in  the  phase  space  plot. 


***  FEL  Phase  Space  Evolution  .*** 
?-*■  ao=3.14;  vQ=0  N=100 
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Figure  2.  Phase-Space  plot  for  an  electron  beam  injected  at  resonance. 

This  clockwise  rotation  can  be  understood  by  looking  at  the  pendulum  equation, 
Eq.  2.28,  and  comparing  it  to  a  mechanical  pendulum.  Each  electron  is  injected  into 
the  undulator  with  an  initial  position  in  phase  space,  (£0,v0).  The  points  (-rc/2,0)  and 
(3tc/2,0)  are  unstable  fixed  points.  This  is  similar  to  the  mechanical  pendulum  being  at 
the  top  of  its  arc.  The  position  (n/ 2,0)  is  a  stable  fixed  point  and  it  represents  the 
mechanical  pendulum  being  at  the  bottom  of  it’s  arc.  The  separatrix,  represented  on 
the  phase  space  plot  by  a  solid  line,  shows  the  boundary  between  open  and  closed 
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orbits  for  the  electrons.  The  open  orbits  represent  a  mechanical  pendulum  with 
enough  total  energy  to  rotate  through  the  unstable  fixed  points.  The  equation  for  the 
separatrix  is  given  by  [4] 

vs2  =  2  la  1  [1-sin(Cs  +<j>)]  ,  (2.43) 

which  results  in  a  total  height  of  4lal1/2.  From  this,  two  important  facts  with  regards  to 
closed  orbits  can  be  seen.  First,  as  the  optical  field  strength  grows,  the  maximum 
height  of  the  separatrix  grows.  Also,  since  the  equation  for  the  separatrix  depends  on 
the  combined  electron  phase,  £  and  the  optical  phase  <j>,  the  separatrix  will  shift  in  £  as 
the  phase  of  the  optical  wave  evolves  in  the  undulator.  This  will  typically  happen  with 
large  current  densities,  j  >n  and  strong  optical  fields,  la  I  >ji. 

In  Fig.  2,  the  electron  beam  was  considered  monoenergetic,  and  thus,  all  of  the 
electrons  are  injected  into  the  undulator  with  the  same  phase  velocity,  v0=0. 
However,  since  the  electron  beam  has  a  finite  length  much  greater  than  the  optical 
wavelength  with  a  relatively  large  number  density  of  electrons  in  the  beam,  it  is 
considered  that  the  electrons  are  injected  with  random  phases  with  respect  to  the 
optical  wave.  In  the  phase  space  plots  used,  the  phases  were  equally  distributed 
between  -n/2  to  3n/2.  From  this  and  Fig.  2,  it  is  seen  that  the  same  number  of 
electrons  will  lose  energy  as  gain,  thus  the  overall  energy  transfer  from  the  electrons 
to  the  optical  field  is  zero. 

ln_Fig.  2,  to  the  right  of  the  phase-space  plot  are  two  plots.  The  top  one  is  a  plot 
of  the  gain  experienced  by  the  optical  field  as  it  propagates  down  the  undulator.  The 
bottom  one  shows  the  evolution  of  the  phase  of  the  optical  wave  as  it  interacts  with 
the  electron  beam  in  the  undulator.  The  gain  is  defined  as  a  function  of  x  and  v0  by 

G  (v0 ,  x)  s  ■—  ,  (2.44) 

ao 

representing  the  fractional  increase  in  the  power  contained  in  the  optical  field,  where 
a0  is  the  initial  optical  field  at  x  =  0.  For  low  current  FELs,  the  average  energy  lost  by 
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one  electron  to  the  optical  field  is  ymc2(<v>- v0)/4tt;/V)  where  <v>  is  the  average 
phase  velocity  of  the  electrons  at  x.  In  weak  optical  fields  the  gain  is  [4] 


<3(v0,t)  =  j 


2-2cos(v0T)-v0TSin(vo  x) 


(2.45) 


where  v0  is  the  initial  phase  velocity  of  the  incident  electron  beam. 


The  bottom  plot  on  the  right-hand  side  of  Fig.  2  represents  the  evolution  of  the 
optical  field  phase,  <j>,  as  a  function  of  x.  From  Eq.  2.23  y/y  ~  0  for  electrons  near 
C +<i>  =  rc/2,and  there  is  little  energy  transferred  between  the  optical  field  and  the 
electrons.  However,  the  interaction  between  the  electrons  and  the  optical  field  results 
in  the  phase  of  the  optical  field  changing  as  shown  in  Equations  2.40  and  2.41 .  This 
optical  phase  shift  affects  the  phase  relationship  between  the  electrons  and  the  optical 
field  as  will  be  shown  in  the  case  of  strong  optical  fields. 

Figure  3  is  a  graph  of  G(v0)  at  x=  1  for  low  current,  j  =  1,  and  weak  field, 
a0  <  7i,  conditions  with  perfect  beam  injection.  This  clearly  shows  that  as  the  initial 
phase  velocity  increases  above  resonance,  more  energy  is 
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Figure  3.  Gain  spectrum  for  a  low  current,  weak  field  FEL. 


transferred  to  the  optical  field.  However,  if  the  initial  phase  velocity  is  too  high,  then  a 
larger  fraction  of  the  electrons  are  no  longer  in  closed  orbits  and  will  not  couple  with 
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the  optical  field.  This  results  in  a  degradation  in  the  gain  as  v0  increases.  Maximizing 
G(v0,t)  with  respect  to  v0  results  in  G  =  0.135y  at  v0  =  2.6.  Note  that  since  G(v0,x) 
is  proportional  to  j,  the  terms  low  gain,  and  low  current  are  interchangeable. 

In  Fig.  4,  the  electron  beam  is  injected  into  the  undulator  off-resonance  at 
v0  =  2.6  where  the  gain  is  a  maximum.  This  corresponds  to  a  higher  average  energy 
in  the  electron  beam  than  in  Fig.  2,  and  more  energy  will  be  transferred  to  the  optical 
field  by  the  electrons  than  will  be  absorbed.  Note  that  some  of  the  electrons  near  the 
edges  of  the  phase  space  plot  will  have  open  orbit  trajectories,  while  those  electrons 
within  the  separatrix  will  have  closed  orbits.  The  electrons  in  closed  orbits  begin  to 
"bunch"  near  £  =  n.  From  the  wave  equation,  Eq.  2.41 ,  it  can  be  seen  that  bunching 
at  this  phase  will  drive  the  optical  field  amplitude  and  produce  gain.  As  the  number  of 
electrons  with  the  same  relative  phase  increases,  the  radiation  from  the  electrons 
becomes  coherent. 


-n/2  £  3n/2  Ox  1 


Figure  4.  Phase-Space  plot  for  an  electron  beam  injected  at  v0  =  2.6. 
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In  order  to  understand  the  interaction  between  the  electron  beam  and  the  optical 
field,  the  pendulum  equation  must  first  be  revisited.  Using  the  definitions  v  =  \  and 
C  =  {k+k0)z -at,  and  assuming  weak  fields  and  low  gain,  the  pendulum  equation  can 
be  expanded  to  get  an  electron’s  phase  and  phase  velocity  with  respect  to  the  optical 
wave  as  a  function  of  x, 


CM  =  Co+voT - |[cos(C0+v0x)-cos(Co+v0x)sin(Co)]  + . 

vo 


(2.46) 


and 


v(x)  =  v0+-^jsin(Co+v0x)-sin(C0) 


ao 


(2.47) 


4  [cos(2^0  +  2v0x)  - cos(2£0)] + cos(v0x)  - 1  - v0xsin(£0)cos(£0  +v0~) 


Using  Eq.  2.46  and  Eq.  2.47,  Eq.  2.41  can  be  solved  for  the  optical  field  strength  and 
phase  evolution  as  a  function  of  x, 


la(x)l  =  a0 


i+y 


2  -  2cos(v0x)  -  v0xsin(v0x) 


2v03 


(2.48) 


and 


H-)  =  j 


2sin(v0x)-v0x(1  +cos(v0x)) 


2Vn 


(2.49) 


For  electrons  injected  with  a  phase  velocity  v0  =  n,  as  in  Fig.  4,  Eq.  2.48  shows  that 
the  optical  field  is  amplified  at  x  =  1 .  With  electrons  injected  at  v0=  -n,  the  optical  field 
is  not  amplified  and,  in  fact,  begins  to  lose  energy  to  the  electrons  as  they  bunch  near 
C  =  0.  When  the  electrons  are  injected  on  resonance,  v0  =  0,  at  low  current  and  weak 
fields,  as  in  Fig.  2,  the  electrons  bunch  near  £  =  n/2.  This  results  in  very  little  gain  in 
the  optical  field  as  shown  by  Eq.  2.48,  however,  Eq.  2.49  shows  that  the  optical  phase 
is  advanced.  With  this  information,  and  Fig.  4,  the  range  of  vq  for  good  coupling 
between  the  electron  beam  and  the  optical  field  is  found  to  be  in  the  range  of  lv0l  <  n. 
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The  FEL  natural  gain  bandwidth  can  be  found  from  Av0  «  4rcA/Ay/y  and  is  given  by 


Ay  „  1 
y  2  N  ' 


(2.50) 


H.  HIGH  CURRENT,  HIGH  GAIN  FELS 

In  the  previous  discussion,  low  currents,  j<  1,  and  weak  fields,  la0i<  n,  were 
assumed  in  order  to  derive  analytical  approximations  which  describe  the  evolution  of 
the  electrons  and  the  optical  field  as  a  function  of  t.  However,  if  j»  1,  these 
assumptions  are  no  longer  valid  and  in  order  to  study  the  FEL  in  these  conditions,  the 
pendulum  and  wave  equations  must  be  used  as  a  starting  point.  The  basic  principle 
of  the  FEL  does  not  change.  The  gain  still  depends  on  the  ability  of  the  FEL  to  bunch 
electrons  and  produce  coherent  radiation. 

Figure  5  is  a  phase  space  plot  of  an  FEL  with  a  dimensionless  current  density  of 
j  =  100,  a  relatively  weak  field  of  a0  =n,  no  spread  in  phase  velocity,  and  uniform 
phase  distribution.  The  initial  phase  velocity  of  vq  =  1 .8  is  the  phase  velocity  at  which 
maximum  gain  is  achieved  for  this  FEL.  At  the  beginning  of  the  undulator,  near  t  =  0, 
the  electrons  begin  to  bunch  near  £+<|>  =  tc/2.  During  this  initial  interaction  time,  there 
is  no  gain  in  the  optical  field,  or  change  in  the  optical  phase,  as  can  be  seen  by  the 
gain  and  phase  plots  on  the  right  in  Fig.  5.  Since  the  electrons  are  near  resonance, 
they  begin  to  drive  the  optical  phase  instead  of  the  optical  field  which  shifts  the 
separatrix  to  the  left.  As  the  optical  phase  shifts,  the  combined  phase  £+<)>  begins  to 
shift  towards  n,  and  the  optical  field  is  amplified.  At  x  ~  0.5,  the  gain  grows 
exponentially  and  continues  to  the  end  of  the  undulator.  The  strong  shift  in  phase  of 
the  optical  field  shifts  the  separatrix  to  the  left  by  ~  n/2  and  the  electrons  bunch  near 
C  =  n/2,  or  £+<}>  «  n  for  optimum  gain. 

The  analytical  solution  for  the  high  current  case  may  be  derived  from  Eq.  2.28 
and  Eq  2.41  by  using  perturbation  methods  and  integrating  over  all  initial  phases  [4], 
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Figure  5.  Phase-Space  plot  for  a  high  current  FEL. 


In  doing  so,  it  is  found  that  for  small  times,  t  « 1,  a(x)  =  aQ(  1  +  //x3/12  +  •  •  -  )  and 
there  is  little  gain  in  the  optical  field.  The  optical  phase,  however,  increases  slowly 
and  in  proportion  to  t3.  For  x<  (2/y)1/3,  the  electrons  are  bunching.  As  soon  as 
bunching  is  sufficient  to  produce  coherent  stimulation,  the  high  current  causes 
exponential  growth  in  the  optical  field  strength.  For  high  current,  j  »  1,  the  optical 
field  strength,  phase,  and  gain  are  given  by 


I  a  (t)  I  ~  -1°  e(//2)'^3T/2  > 

3 


<Kt)  ~  C/72)1/3x/2  , 


(2.51) 


(2.52) 


G(t)  =  —  e^72)1^ 

9 


(2.53) 


Figure  6  shows  the  high  gain  spectrum  for  j  =  100  and  a0  =  1 .  The  gain 
spectrum  is  no  longer  antisymmetric  about  v0  =  0,  instead  it  is  nearly  symmetric  about 
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the  phase  velocity  which  produces  the  maximum  gain.  This  is  true  in  general  for  high 
current  FELs.  The  peak  gain  occurs  for  a  phase  velocity  that  is  closer  to  resonance 
than  in  the  low  current  case.  Gain  continues  to  be  positive  even  for  negative  values  of 
v0  near  resonance.  The  gain  spectrum  bandwidth  at  which  G(v0)  =  1/e  Gmax,  is 
Av0  =  4/ 1/6  ~  2 n  and  is  much  wider  than  that  for  the  low  current  FEL. 


j=10G  a  =1 

...  o 

71 
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The  characteristic  energy  spread  at  which  gain  degradation  begins  no  longer 
follows  that  of  the  low  current  FEL  of  oG  <  n.  The  FEL  becomes  more  resistant  to 
gain  degradation  when  the  optical  field  amplitude  grows  exponentially.  From  Eq.  2.53 
it  can  be  seen  that  for  high  current  FELs,  the  growth  rate  is  exponential  for  strong 
fields.  -The  characteristic  energy  spread  at  which  gain  degradation  begins  is 
comparable  to  this  growth  rate,  oG  ~  y'1/3  [4]. 

I.  SATURATION  IN  STRONG  OPTICAL  FIELDS 

So  far,  only  weak  optical  fields  have  been  considered.  Equations  2.28  and  2.41 
are  both  valid  for  all  cases  of  strong  and  weak  fields  as  well  as  low  and  high  current 
FELs.  Using  the  phase-space  plots  will  help  in  understanding  the  interactions 
associated  with  strong  optical  fields  and  the  electron  beam.  Figure  7  is  a  phase-space 
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plot  for  a  low  current,  j  =  1 ,  FEL  with  initial  optical  field  strength  of  a0  =  20.  The 
electron  bunching  seen  in  Fig.  4  has  been  extended  until  the  condition  of 
overbunching  is  reached.  The  electrons  began  to  bunch  at  C+<t>  =  n.  As  they  continue 
down  the  undulator,  they  transfer  more  energy  to  the  optical  field  until  they  have 
transferred  the  maximum  amount  of  energy  to  the  optical  field.  At  this  point,  they 
begin  to  absorb  energy  from  the  optical  field.  This  is  called  overbunching  and  is 
indicated  by  a  reduction  in  the  gain.  The  gain  plot  of  Fig.  7  shows  the  initial  rise  in 
gain  due  to  the  strong  optical  field,  and  at  1/2,  maximum  bunching  occurs  as 
indicated  by  the  maximum  rate  of  change  in  gain.  At  t  =  0.8,  the  electrons  have 
transferred  all  of  their  available  energy  to  the  optical  field  and  begin  absorption.  The 
final  gain  of  0.03  is  significantly  reduced  from  that  of  0.135  in  Fig.  4.  This  reduction  in 
gain  in  the  presence  of  strong  optical  fields  is  expected  for  two  reasons.  In  order  to 
achieve  the  same  gain  with  a  larger  initial  optical  field  strength,  the  amount  of  energy 
transferred  to  the  optical  field  must  also  be  proportionally  larger.  The  overbunching  of 
the  electrons  inhibits  the  ability  of  the  electron  beam  to  transfer  energy  to  the  optical 
field.  The  evolution  to  an  overbunching  condition  necessarily  requires  that  the 
electrons  pass  through  the  condition  where  they  are  bunched  at  £+<|>  =  n/2  where  they 
will  advance  the  phase  of  the  optical  field.  This  can  be  seen  in  Fig.  7  where  the 
optical  phase  shift  increases  near  x  »  0.5.  This  particular  case  does  not  result  in  a 
significant  phase  shift  in  the  optical  field,  <j>  <  0.003,  since  the  electron  current  is  low. 

Figure  8  is  a  phase-space  plot  of  an  FEL  with  high  current  and  strong  optical 
fields.  The  overbuncing  is  noticeable.  The  optical  phase  also  advances  more  than  in 
the  low  current  case.  This  shifts  the  separatrix  to  the  left.  The  overall  gain  is,  as 
expected,  significantly  less  than  for  the  weak  field  case  of  Fig.  5. 

The  evolution  to  saturation  in  an  FEL  oscillator  from  weak  fields  is  accompanied 
by  a  change  in  the  optical  frequency.  In  an  oscillator,  the  semi-transparent  mirror  that 
allows  the  laser  light  to  leave  the  oscillator  cavity  results  in  a  reduction  in  the  optical 
field  strength.  This  loss,  combined  with  other  losses  in  the  oscillator  is  quantified  by 
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Figure  7.  Phase-Space  plot  for  a  low  current,  strong  field  FEL. 

the  resonator  Q  factor  defined  in  such  a  way  that  in  the  absence  of  an  electron  beam, 
the  optical  field  is  reduced  by  a  factor  of  e~VQ  each  pass.  If  the  gain  produced  by  the 
electron  beam  is  sufficient  to  make  up  for  these  losses,  then  the  optical  field  will  grow, 
in  weak  fields,  mode  competition  results  in  the  optical  frequency  being  centered 
around  the  frequency  corresponding  to  v0  =  2.6.  As  the  optical  field  grows,  the  phase 
velocity  for  maximum  gain  shifts  to  the  right  and  the  optical  wavelength  decreases.  In 
weak  optical  fields,  with  low  Q,  the  FEL  gain  may  not  be  sufficient  to  overcome  the 
losses  and  the  steady  state  optical  field  is  zero.  For  larger  values  of  Q,  the  FEL  gain 
is  sufficient  to  overcome  the  losses  and  the  optical  field  grows.  When  the  optical  field 
reaches  the  point  where  saturation  occurs  during  each  pass,  then  the  gain  is  reduced 
due  to  this  saturation.  As  the  gain  falls  off,  the  optical  field  continues  to  grow  until  the 
gain  is  insufficient  to  make  up  for  the  losses  in  the  cavity  and  the  growth  rate  of  the 
optical  field  is  reduced.  The  combination  of  saturation,  cavity  losses  and  gain 
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Figure  8.  Phase-Space  plot  for  a  high  current,  strong  field  FEL. 

work  together  to  bring  the  FEL  into  steady  state  operation.  This  typically  takes  several 
hundred  passes.  This  concept  is  important  in  the  study  of  SELENE. 

J.  OPTICAL  KLYSTRON  FELS 

One  of  the  beneficial  features  of  the  FEL  is  the  ability  to  change  the 
characteristic  output  of  the  FEL  by  changing  the  undulator  design.  Several  undulator 
designs  are  used  for  altering  the  interaction  between  the  electron  beam  and  the  optical 
field.  One  of  these  designs  shown  in  Fig.  9  is  the  optical  klystron.  The  optical  klystron 
is  used  to  provide  increased  gain  in  weak  fields.  This  is  accompanied  by  saturation  at 
lower  optical  field  strengths  compared  to  the  standard  FEL.  The  optical  klystron 
consists  of  an  FEL  oscillator  constructed  with  a  dispersive  or  drift  section,  positioned  at 
x  =  0.5,  between  two  undulator  sections  . 

The  first  undulator  section,  or  modulator,  is  used  to  prepare  the  electrons  for 
bunching.  When  the  electrons  enter  the  dispersive  section,  the  interaction  with  the 


***  fel  Phase  Space  Evolution  *** 
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Figure  9.  Optical  klystron  configurations. 


optical  field  stops.  In  the  drift  space  ,  this  happens  because  the  electrons  are  far  from 
resonance,  while  in  the  dispersive  section  the  electrons  are  bent  out  of  the  optical  field 
by  a  large  magnetic  field.  In  the  dispersive  section,  the  electrons  are  bunched  due  to 
the  differences  in  phase  velocities  between  the  electrons.  The  electrons  with  phase 
velocities  greater  than  the  average  experience  a  larger  Lorentz  force  than  the  slower 
electrons.  As  they  pass  through  the  magnetic  field,  they  have  a  longer  path  length  to 
travel  before  reentering  the  second  undulator  section.  The  opposite  is  true  for  the 
slower  electrons.  Thus,  as  the  electrons  enter  into  the  second  undulator  section,  or 
radiator,  the  electrons  are  bunched  near  the  same  phase  relative  to  the  optical  field. 
With  the  electrons  bunched  prior  to  entering  the  second  undulator  segment,  the  ability 
of  the  electrons  to  radiate  and  transfer  energy  to  the  optical  field  is  much  greater. 

Since  the  electrons  do  not  interact  with  the  optical  field  while  they  are  in  the 
dispersive  section,  they  do  not  transfer  energy  to  or  from  the  optical  field.  However, 


they  do  experience  a  phase  shift  that  is  directly  proportional  to  their  phase  velocity  and 
the  dispersive  strength,  A C  =  vD.  D  is  the  dispersive  strength  of  the  klystron  given  in 
equivalent  undulator  lengths,  D  =  Nd/N,  where  Nd  is  the  number  of  undulator  periods 
that  would  produce  the  same  bunching.  The  dispersive  strength  is  governed  by  the 
strength  of  the  magnetic  fields  in  the  dispersive  section  and  the  undulator  parameter 
K.  A  dispersive  strength  of  one  is  equivalent  to  having  doubled  the  length  of  the 
undulator.  Equations  2.28  and  2.41  are  still  valid  throughout  the  length  of  the 
undulator  for  an  optical  klystron.  For  jD  « 1  and  in  the  limit  of  D->°°  and  v0  0  with 
the  product  Dv0  fixed,  the  gain  as  a  function  of  v0  is  [7] 

G(v0)  ~  -^sin(voD) ,  (2.54) 

with  a  peak  gain  of  /D/4  at  v0  =  n/2D .  This  is  illustrated  in  Fig.  1 0. 
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Figure  1 0.  Gain  spectrum  for  a  low  current,  weak  field  optical  klystron  FEL. 

Figure  10  shows  that  there  are  several  peaks  in  the  region  lv0l<2n.  From  Eq. 
2.54,  it  is  apparent  that  the  electron  beam  quality  should  be  such  that  Av0<nlD  in 
order  to  maximize  coupling  and  gain.  Compared  to  the  non-klystron  FEL,  where 
Av0  <7i  is  required,  this  is  a  significant  restriction  on  using  an  optical  klystron. 
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Figure  1 1  is  a  phase-space  plot  for  an  optical  klystron  in  weak  fields  and  low 
current.  The  electrons  are  injected  into  the  undulator  with  uniform  phase  velocity  of 
v0  =  1  -1  for  maximum  gain  and  are  represented  in  light  gray.  As  the  electrons 
propagate  down  the  undulator,  they  become  darker  until  they  reach  the  end  of  the 
undulator  where  they  are  shown  in  black.  The  evolution  of  the  electrons  in  the 
modulator  is  similar  to  that  of  Fig.  4.  The  gain  plot  to  the  right  of  the  phase-space  plot 
shows  that  the  optical  field  is  not  driven  until  after  the  dispersive  section  at  t  =  0.5.  At 
this  point,  the  electrons  are  given  a  step  jump  in  phase,  A C  =  v£>,  which  can  be  seen 
in  Fig.  1 1 .  This  results  in  all  of  the  electrons  moving  to  the  right  a  distance  which  is 
proportional  to  the  phase  velocity  of  the  individual  electrons.  The  electron  beam  and 
optical  field  interaction  continues  from  this  point  to  achieve  the  FEL  output.  The 
bunching  is  easy  to  see  at  £ + <j>  =  n  at  the  end  of  the  second  undulator.  The  final  gain 
of  0.04  is  almost  three  times  that  of  the  non-klystron  undulator.  This  shows  the 
usefulness  of  using  the  optical  klystron  configuration  in  FELs  that  have  insufficient  gain 
to  startup. 

In  the  presence  of  strong  optical  fields,  the  characteristics  of  the  optical  klystron 
changes  dramatically.  Figure  1 2  shows  the  gain  spectrum  for  an  optical  klystron  with 
jD  <  1  and  strong  optical  fields.  Figure  13  is  the  phase-space  plot  for  v0  =  4.1  at  which 
the  maximum  gain  is  achieved.  The  positions  of  the  electrons  in  phase-space  at 
t  =  0.5,  just  prior  to  the  dispersive  section  is  plotted  in  light  gray.  The  evidence  of  the 
strong  optical  fields  is  evident  in  that  the  electrons  are  already  bunched  prior  to 
entering  the  dispersive  section.  When  the  electrons  enter  the  dispersive  section,  the 
electrons  are  bunched  at  £+<(>  =  n  and  are  ail  given  a  step  displacement  in  £.  Since 
they  were  previously  bunched,  this  results  in  the  electrons  being  dispersed  instead  of 
being  bunched.  This  results  in  a  reduction  in  gain  by  a  factor  of  ~  7.5  as  compared  to 
the  weak  field  case  of  Fig.  1 1 .  The  electrons  continue  to  interact  with  the  optical  field 
as  they  propagate  down  the  radiator  section  and  the  final  phase-space  positions  are 
plotted  in  dark  gray  in  Fig.  13. 
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Figure  1 1 .  Phase-Space  plot  for  a  low  current,  weak  field  optical  klystron. 
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Figure  12.  Gain  spectrum  for  a  low  current,  strong  field  optical  klystron  FEL. 

In  stronger  optical  fields,  the  modulator  is  efficient  at  bunching  the  electrons,  so 
that  when  the  electrons  reach  the  dispersive  section,  the  AC  induced  by  the  dispersive 
section  essentially  redistributes  the  electrons  evenly  in  phase  with  a  large  energy 
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Figure  13.  Phase-Space  plot  for  a  low  current,  strong  field  optical  klystron  FEL, 
So  =  10- 


spread.  This  results  in  the  electrons  entering  the  radiator  as  a  low  quality  electron 
beam  instead  of  the  highly  bunched  electron  beam  the  dispersive  section  is  designed 
to  produce  in  weak  fields. 

The  result  of  the  dispersion  effects  of  the  optical  klystron  in  strong  optical  fields  is 
evident  in  Fig.  14,  where,  at  x  =  0.5,  the  gain  flattens  out  until  the  electrons  can  be 
rebunched  by  the  radiator  section.  The  final  gain  from  this  FEL  is  significantly  reduced 
to  about  half  that  of  Fig.  13.  This  reduction  in  gain  in  the  presence  of  strong  optical 
fields  is  the  reason  why  dispersive  sections  are  designed  to  be  used  to  produce  high 
gain  only  in  weak  field  FELs.  This  property  of  reduced  gain  in  strong  field  optical 
klystron  FELs  can  actually  become  an  advantage,  and  will  be  exploited  in  the  design 
of  SELENE. 
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Figure  14.  Phase-Space  plot  for  a  low  current,  strong  field  optical  klystron  FEL,  a0  =  20 
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III.  SELENE  AND  NOVOSIBIRSK  PROPOSALS 


A.  INTRODUCTION 

As  a  joint  venture,  Russia  and  the  United  States  are  pursuing  the  development  of 
a  high  power  laser  system  to  beam  energy  to  satellites.  The  Budker  Institute  at 
Novosibirsk,  Russia,  is  building  a  200  kW  system.  The  proposal  for  SELENE  at  China 
Lake,  California,  is  intended  to  expand  on  this  design  to  achieve  a  10  MW  system. 
Figure  15  is  a  schematic  of  the  proposed  system  for  SELENE.  This  consists  of  a 
three-section  optical  klystron  oscillator  followed  by  a  radiator.  The  three  undulator 
sections  and  the  undulator  in  the  radiator  are  identical.  The  electron  beam  will  be 
provided  by  a  51  MeV  racetrack  microtron  recuperator  (RTMR)  [9],  The  oscillator  is 
designed  to  bunch  the  electron  beam  prior  to  entering  the  radiator.  The  dispersive 
sections  in  the  oscillator  are  designed  to  be  variable  in  strength.  The  electron  beam  is 
bent  out  of  the  oscillator  and  sent  to  the  radiator  which  uses  coherent  spontaneous 
emission  from  the  bunched  electron  beam  to  generate  a  high  power  laser  beam.  This 
beam  is  transmitted  along  a  two  mile  vacuum  tube  to  allow  it  to  expand  by  diffraction 
to  reduce  the  incident  laser  intensity  on  the  mirror.  The  beam  is  reflected  off  of  an 
adaptive  mirror  either  to  the  satellites  directly,  or  to  a  beacon  satellite  which  will  then 
direct  the  laser  to  the  proper  satellites. 

B.  SIMULATION  PARAMETERS 

The  electron  beam  is  20-100  ps  in  duration  with  a  repetition  frequency  of  2  -  45 
MHz.  The  peak  current  is  20  -  1 00  A.  The  emittance  is  estimated  at  0.628  mm-mrad. 
The  51  MeV  electrons  have  a  Lorentz  factor  of  y=100.  The  peak  electron  beam 
current  of  100  A,  and  the  beam  area  of  nrb  =  0.002cm2,  where  rb  is  the  electron  beam 
radius,  results  in  an  average  filling  factor  of  F  =  0.18  over  the  entire  system  and 
F  =  0.23  in  the  oscillator.  The  total  dimensionless  current  density  is  the  product  of  the 
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Figure  15.  Schematic  diagram  of  the  SELENE  FEL  system. 


dimensionless  current  density  j  and  the  filling  factor  P.  For  further  ease,  the  total 
dimensionless  current  density  will  be  represented  by  jF  j.  The  parameters  above 
result  in  a  total  dimensionless  current  of  y  =  180  for  the  entire  system,  and  j  =  100  for 
just  the  oscillator.  The  design  energy  spread  is  Ay/y  =  0.045%  which  results  in  an 
initial  phase  velocity  spread  of  Av  =  aG  =  4nA/Ay/y  =  0.67  for  the  oscillator  and 
aG  =  0.9  for  the  system  as  a  whole.  Assuming  the  electron  beam  has  been  groomed 
for  matched  angular  and  radial  distributions,  the  rms  angular  deviation  from  the  z  axis 
is  0  =  ejnrb  =  0.00025  and  the  angular  spread  is  characterized  by  o0  =  0.32. 

Each  undulator  section  consists  of  40  undulator  periods.  Each  period  is  9  cm  in 
length,  for  a  total  length  of  3.6  m  for  each  section.  The  proposed  dispersive  section 
strength  is  equivalent  to  D  =  0.5  for  each  dispersive  section  in  the  oscillator.  The 


40 


space  between  the  oscillator  and  the  radiator  also  acts  like  a  dispersive  section  with 
an  equivalent  dispersive  strength  of  Dr  =  0.25.  When  the  oscillator  is  studied  as  a 
separate  unit,  the  dispersive  strength  for  each  dispersive  section  scales  to  D  =  0.67. 
During  the  study  of  the  oscillator  it  is  desirable  to  compare  the  effects  of  this  three- 
section  oscillator  to  a  two-section  oscillator.  For  the  three-section  oscillator,  the 
dispersive  strength  stated  is  for  one  of  the  dispersive  sections.  For  the  two-section 
oscillator  used  to  compare  with  SELENE’s  oscillator,  the  dispersive  strength  of  D  =  1 .3 
will  be  used  to  closely  approximate  the  total  dispersive  strength  of  both  dispersive 
sections  in  the  three-section  klystron.  The  mirrors  in  the  oscillator  section  are 
separated  by  a  distance  of  79  m.  The  FEL  undulator  parameter  of  K  =  1.4  [9] 
corresponds  to  a  optical  wavelength  of  13  p.m  and  £  =  0.3,  where  £,  is  defined  following 
EQ.  2.30. 

C.  OPTICAL  DEVICE  LIMITATIONS 

A  major  concern  in  the  SELENE  FEL  design  is  the  durability  of  the  mirrors  in  the 
oscillator.  The  beam  from  the  radiator  has  a  two  mile  propagation  distance  to  spread 
out  and  reduce  the  laser  intensity  on  the  mirrors  that  project  the  laser  to  the  satellites. 
The  oscillator,  on  the  other  hand,  has  an  11  meter  undulator  section  in  a  79  meter 
oscillator  cavity.  The  high  electron  beam  current  results  in  high  gain  and  thus  strong 
optical  fields  after  saturation  has  occurred.  Mirror  construction  is  a  highly  specialized 
field  and  the  materials  used  to  construct  them  vary  widely  in  their  ability  to  withstand 
high  intensity  laser  light.  Mirrors  are  available  to  operate  with  peak  intensities  of 
1  GW /cm 2  with  pulse  lengths  of  0.1  ns  and  a  ten  percent  duty  cycle  [10].  Mirrors  have 
been  demonstrated  which  can  withstand  an  average  power  density  of  10-100kW/cm2 
with  a  continuous  duty  cycle  [11].  In  order  to  provide  a  safety  margin,  the  maximum 
intensity  allowed  on  the  mirrors  is  assumed  to  be  50kW/cm2  =  5X1 011  ergs /cm  2s  for 
SELENE.  If  the  Rayleigh  length  is  matched  to  an  individual  undulator  section,  then 
z0~  4  meters.  The  expansion  of  the  optical  wave  due  to  diffraction  is  given  by 
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(3.2) 


[  C2 
=  1+% 
l  Z0. 

where  Sm  is  the  mirror  separation  distance,  w  is  the  radius  of  the  beam  spot  size  at 
the  mirrors,  w0  is  the  optical  beam  radius  at  the  waist  which  is  usually  in  the  center  of 

the  undulator,  and  z0X  =  nw2°.  The  magnitude  of  the  optical  field  intensity  at  the 
mirrors  is  [12] 

I5>l  =  l£l2c/4rc  .  (3.3) 

Thus  the  maximum  allowable  optical  field  strength  in  the  undulator  to  prevent  mirror 
damage  is 


=  1+- 


4nNeK(J  0(5)~ 
y2 me 2 


Using  the  numbers  presented  above  for  this  system,  the  maximum  dimensionless 
optical  field  strength  in  the  undulator  is  la  I  »  30.  This  is  well  into  the  strong  field 
regime.  Once  the  optical  field  is  reduced  below  this  value  for  steady  state  operation, 
the  concern  is  to  ensure  the  oscillator  is  able  to  perform  the  function  for  which  it  was 
designed,  that  is,  to  provide  a  bunched  electron  beam  to  the  radiator.  As  discussed 
previously,  strong  optical  fields  result  in  overbunching  and  thus  this  will  inhibit  the 
ability  of  the  radiator  to  produce  coherent  spontaneous  emission.  The  spontaneous 
emissions  are  required  to  establish  an  optical  field  for  the  electron  beam  to  interact 
with  in  the  single  pass  radiator. 


Figures  13  and  14  demonstrate  that  strong  optical  fields  are  undesirable  if  the 
electrons  are  required  to  be  bunched  prior  to  entering  the  SELENE  radiator.  In  strong 
optical  fields,  the  electrons  are  bunched  prior  to  entering  the  first  dispersive  section. 
When  this  happens,  the  dispersive  section  disperses  the  electrons  instead  of  bunching 
them.  It  is  important  to  remember  that  the  electron  beam  for  SELENE  will  encounter  a 
third  dispersive  section  prior  to  entering  the  radiator. 
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D.  DISPERSIVE  SECTION  EFFECTS 

The  oscillator  design  of  the  SELENE  system  provides  some  challenging  goals  in 
trying  to  reduce  the  optical  field  strength  sufficiently  to  prevent  damage  to  the  mirrors 
and  to  prevent  destroying  the  bunching  of  the  electrons  prior  to  entering  the  radiator. 
Figure  16  shows  the  result  of  using  optical  klystrons  in  strong  optical  fields.  This 
figure  was  produced  using  a  series  of  gain  spectra  similar  to  Fig.  1 2  where  the  initial 
optical  field  strength  was  varied  and  the  peak  gain  was  recorded.  In  the  weak  field 
regime,  a0  <  n,  the  dispersive  sections  increase  the  gain  significantly.  This  will  result 

in  the  optical  field  growing  rapidly  during  each  pass  until  the  strong  field  condition  is 
reached. 


Initial  Dimensionless  Optical  Field  Strength 


Figure  16.  Single  pass  maximum  gain  vs.  dimensionless  optical  field  strength  for  a 
three-section  klystron,  a  two-section  klystron,  and  a  single  undulator  FEL 
oscillator. 

Even  though  the  two  dispersive  sections  in  the  three-section  klystron  have  the 
same  total  dispersive  strength  as  the  single  dispersive  section  in  the  two-section 
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klystron,  the  three-section  klystron  produces  more  gain  in  weak  fields  than  the  two- 
section  klystron.  However,  as  the  optical  field  strength  increases,  the  gain  in  the 
three-section  klystron  oscillator  becomes  less  than  the  gain  for  the  two-section 
klystron.  As  the  optical  field  continues  to  grow,  the  gain  for  the  three-section  klystron 
produces  is  reduced  below  that  of  the  single  undulator  oscillator.  Figure  16  is  based 
on  the  single  pass  gain  in  an  FEL  with  the  same  operating  parameters  as  SELENE’s 
oscillator. 

Figure  17  is  a  plot  of  the  single  pass  maximum  gain  as  a  function  of  the 
dispersive  strength  with  an  initial  optical  field  strength  of  a0  =  10.  Three  important 
features  are  visible  in  this  plot.  First,  for  both  the  three  section  and  the  two  section 
klystrons,  increasing  the  dispersive  strength  beyond  D  =  2  has  little  effect  on  the 
maximum  gain  of  the  system  for  D  <,4.  Next,  it  is  evident  that  the  three-section 
klystron  has  less  gain  at  a0  =  1 0  than  the  two  section  klystron  for  D  £  0.3.  This 
agrees  with  Fig.  16.  The  other  major  feature  of  note  is  the  peak  in  gain  at  D  «  0.5. 
The  origin  of  Fig.  17  is  a  series  of  gain  spectrum  plots  similar  to  Fig.  12.  The 
maximum  gain  for  each  dispersive  strength  is  plotted  in  Fig.  17.  At  D  s  0.5  for  both 
systems,  the  initial  phase  velocity  giving  the  maximum  gain  shifts  to  a  higher  phase 
velocity. 

Figure  18  shows  this  phenomenon  in  a  form  easier  to  interpret.  If  the  phase 
velocity  is  kept  constant  and  <jq  is  small  enough  to  ensure  that  few  electrons  are 
outside  the  central  maximum  shown  in  Fig.  18,  then  the  gain  will  drop  rapidly  as  the 
dispersive  strength  is  increased.  However,  if  the  electron  phase  velocity  is  maintained 
at  that  velocity  for  which  maximum  gain  is  achieved,  then  the  result  is  the  effect  shown 
in  Fig.  17.  The  process  of  natural  mode  competition  in  the  FEL  oscillator  maintains 
the  relative  phase  velocities  between  the  optical  field  and  the  electron  beam  so  that 
the  average  phase  velocity  of  the  electron  beam  is  that  which  produces  the  maximum 
gain  in  the  optical  field. 
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Figure  17.  Single  pass  maximum  gain  vs.  dispersive  strength  for  a  three-section 
klystron  and  a  two-section  klystron  FEL  with  an  initial  dimensionless 
optical  field  strength  of  a0  =  10. 

So  far,  only  the  effects  of  dispersive  sections  on  the  maximum  gain  in  a  single 
pass  have  been  shown.  In  Fig.  16,  the  cumulative  effects  of  multiple  passes  may  be 
extrapolated.  As  the  optical  field  strength  increases,  each  of  the  three  example 
oscillators  have  a  different  rate  at  which  the  gain  is  reduced.  If  Fig.  16  is  expanded 
so  that  the  point  at  which  each  oscillator  has  a  gain  of  unity  can  be  determined,  the 
results  are  not  as  one  might  expect  at  first  glance.  The  three-section  klystron  is  the 
first  oscillator  to  reach  a  maximum  gain  of  unity  as  the  optical  field  is  increased.  The 
oscillator  without  the  dispersive  section  is  next  to  reach  a  gain  of  unity.  The  slope  of 
the  curve  for  each  oscillator  design  is  different  as  it  crosses  the  G  =  1  line  and  is  not 
constant  as  the  optical  field  increases.  The  region  about  which  the  oscillators  reach 
G  =  1  is  between  a0  =  44  and  a0  =  51 .  This  would  lead  to  the  conclusion  that  the 


45 


ln(l+G) 


16  2 


Figure  18.  Single  pass  gain  spectrum  as  a  function  of  dispersive  strength,  D,  for  a 
two-section  klystron  FEL. 

steady-state  optical  field  strength  after  many  passes  may  be  significantly  different  in 
each  of  the  oscillators. 

Figures  19,  20,  and  21  show  the  evolution  of  the  optical  field  after  2000  electron 
pulses  have  been  passed  through  each  oscillator.  Two  new  FEL  parameters  are 
introduced  at  this  point,  desynchronism  and  slippage  distance.  Desynchronism  and 
slippage  distance  may  be  understood  by  following  the  optical  field  through  a  full  pass 
through  the  undulator  from  the  point  at  which  it  encounters  the  first  electron  pulse  to 
the  next  electron  pulse.  When  the  electrons  have  a  phase  velocity  of  v0  =  0,  they  are 
said  to  be  at  resonance.  This  means  that  for  every  undulator  period,  they  slip  back 
along  the  light  wave  by  one  wavelength  of  light.  This  results  in  the  electron  beam 
falling  behind  the  optical  field  by  a  distance  of  NX  as  it  travels  the  entire  length  of  the 
undulator.  The  normalized  length  of  the  electron  pulse  is  oz  =  le/NX,  where  le  is  the 
length  of  the  electron  pulse.  Thus  a  slippage  distance  characterized  by  oz  <  1  means 
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that  the  electron  beam  has  slipped  back  along  the  optical  field  a  distance  that  is 
greater  than  the  length  of  the  electron  micropulse.  This  results  in  distortion  in  the 
shape  of  the  optical  pulse  since  the  electron  beam  is  not  interacting  with  the  same 
point  in  the  optical  field  as  it  travels  down  the  undulator.  As  the  electron  beam  slips 
back  along  the  optical  field,  the  optical  field  grows  stronger  farther  back  along  the 
length  of  the  optical  pulse.  If  a  correction  is  not  implemented,  this  will  result  in  an 
optical  pulse  skewed  towards  the  beginning  of  the  undulator.  This  effect  is  corrected 
by  introducing  desynchronism.  A  desynchronism  of  d  =  1  would  mean  that  the 
electron  pulse  is  injected  into  the  undulator  so  that  it  is  one  slippage  distance  ahead  of 
the  optical  pulse  relative  to  the  position  it  was  introduced  into  the  previous  pass. 
Positive  desynchronism  is  accomplished  by  increasing  the  length  of  the  resonator 
cavity  by  a  small  fraction  of  the  slippage  distance  so  that  the  optical  pulse  arrives  late 
with  respect  to  the  electron  micropulse.  Values  of  desynchronism  are  typically  less 
than  1%  of  the  slippage  distance.  Negative  values  of  desynchronism  are  used  in 
some  applications  [1 0]. 

Figures  19,  20  and  21  contain  a  large  amount  of  information  in  ten  different  plots. 
The  plot  on  the  bottom  left  gives  an  indication  of  how  far  the  electron  beam  slips  back 
along  the  optical  field  in  each  pass  relative  to  the  length  of  the  electron  pulse  where  z 
is  measured  relative  to  the  optical  field  in  slippage  distances.  The  dark  curve  is  the 
electron  pulse  at  the  beginning  of  the  undulator  and  the  light  gray  shape  is  the 
electron- pulse  at  the  end  of  the  undulator.  The  difference  in  the  position  of  the 
electron  pulse  from  the  beginning  to  the  end  of  the  undulator  is  Az  =  1 ,  the  slippage 
distance,  or  approximately  1/az  ~  25%  of  the  length  of  the  electron  pulse  for  SELENE. 
The  plot  in  the  center  on  the  bottom  is  the  weak  field  gain  spectrum  for  the  initial 
conditions  of  the  oscillator.  The  plot  on  the  bottom  right  is  the  power,  P,  in  the  optical 
field  at  the  end  of  each  pass.  The  power  is  the  average  of  la(z)l2  over  the  entire 
window  shown  in  the  upper-left  corner.  When  this  plot  is  level  over  many  passes,  the 
oscillator  may  be  considered  to  be  in  steady  state. 
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*  SELENE  FEL  Pulse  Evolution  * 
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Figure  19.  FEL  output  parameters  after  2000  electron  pulses  for  a  three-section 
optical  klystron  oscillator. 

The  three  middle  plots  show  the  development  of  the  FEL  over  many  passes. 
The  left  plot  demonstrates  how  the  optical  field  develops.  The  striated  features 
indicate  that  the  shape  of  the  optical  field  is  no  longer  smooth.  The  middle  plot  shows 
the  power  distribution  as  a  function  of  the  phase  velocity  of  the  electrons  and  the 
number  of  passes.  The  right  graph  shows  the  distribution  of  the  electrons  at  the  end 
of  each  pass  as  a  function  of  phase  velocity.  From  this,  it  can  be  seen  that  the 
oscillator  with  the  single  dispersive  section  is  more  efficient  than  the  other  two 
configurations  at  bunching  the  electrons. 

The  three  plots  at  the  top  of  Figs.  19,  20  and  21  show  the  optical  field,  power, 
and  electron  distribution  at  the  end  of  the  last  pass.  The  two  symbols  at  the  top  of  the 
power  and  electron  distribution  graphs  show  the  average  distribution  of  these  two 
functions.  The  triangle  indicates  the  average  for  the  first  electron  pulse.  The  line 
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Rgure  20.  FEL  output  parameters  after  2000  electron  pulses  for  a  single  undulator 
oscillator. 

indicates  the  average  at  the  end  of  the  last  pass.  The  scale  figure  for  the  optical  field 
strength,  located  in  the  upper-left  corner,  indicates  that  the  maximum  optical  field 
strength  is  plotted  in  light  gray,  while  the  minimum  optical  field  is  plotted  in  black.  This 
scale  applies  for  the  power  and  electron  distribution  graphs,  but  it  is  normalized  for 
each  application. 

Figures  19,  20  and  21  indicate  that  if  the  damage  threshold  of  the  mirrors  is 
sufficient,  the  SELENE  parameters  could  provide  a  peak  dimensionless  optical  field 
strength  greater  than  700.  This  results  in  a  peak  power  density  in  excess  of 
24  MW /cm2.  The  three-section  klystron  does  have  the  lowest  peak  optical  field 
strength,  while  the  undulator  without  the  dispersive  section  has  the  highest.  The  value 
of  the  optical  field  given  by  the  scale  is  the  maximum  optical  field  strength  developed 
during  the  entire  run  at  any  position  along  the  optical  wave.  It  is  very  easy  to  see  that 


49 


11? ;;  SEIBNtE  F§£  P&Lse 
j=I0a:irCTz*3.i-  ae=0.32  0rG=tGi67  d=Q.;<30I 
Q-20  D*1.3  N~120 


-4  z  4-101  v  101-101  v  i01 


-4  z  4-101  v  1010  n  2000 


Figure  21.  FEL  output  parameters  after  2000  electron  pulses  for  a  two-section 
optical  klystron  oscillator. 

the  shape  of  the  optical  wave  is  definitely  not  a  uniform  or  smooth  waveform.  This  is 
due  in  particular  to  the  high  optical  field  strength  and  large  electron  current  interacting 
in  the  undulator.  The  process  of  over  bunching  results  in  small  pockets  of  electrons 
along  the  length  of  the  optical  wave  which  produce  local  amplifications  much  greater 
than  the-  average,  and  thus  the  sharp  peaks  are  found  along  the  optical  wave.  The 

presence  of  the  dispersive  section  in  strong  optical  fields  tends  to  aggravate  this 
problem. 

E.  SELENE  OPTIMIZATION 

From  Figs.  19,  20,  and  21,  it  is  easily  seen  that  the  current  design  of  the 
SELENE  project  produces  too  much  gain  in  the  oscillator,  as  the  optical  field 
produced  at  steady  state  is  much  too  high  (a  «  750)  for  current  mirror  technology  as 
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well  as  producing  an  ill-conditioned  electron  beam  for  the  radiator.  Guided  by  the 
basic  principles  of  FEL  design,  it  is  desired  to  retain  the  basic  SELENE  structure  while 
finding  ways  to  reduce  the  optical  field  in  the  oscillator.  This  must  be  accomplished 
while  maintaining  the  function  of  the  oscillator,  which  is  to  prepare  the  electron  beam 
for  the  radiator.  Several  parameters  may  be  varied  in  the  design  of  the  oscillator 
without  incurring  unusual  deviations  from  the  original  design.  Three  of  these 
parameters  will  be  used  to  find  a  method  which  will  reduce  the  peak  optical  field  in  the 
SELENE  oscillator. 

Desynchronism  is  one  useful  way  of  reducing  the  output  power  of  the  oscillator. 
The  purpose  behind  desynchronism  is  to  try  and  maintain  the  shape  of  the  optical  field 
as  it  propagates  down  the  undulator.  The  top  left-hand  plot  in  Fig.  22  shows  the 
effects  of  electron  slippage  on  the  shape  of  the  optical  pulse.  As  the  electrons  slip 
back  along  the  optical  pulse  each  pass,  the  trailing  edge  of  the  pulse  is  amplified  more 
than  the  forward  edge. 

Introducing  desynchronism  also  produces  a  reduction  in  gain.  As  the  electrons 
are  moved  ahead  of  the  optical  field  (positive  desynchronism)  or  behind  the  optical 
field  (negative  desynchronism)  some  of  the  electrons  are  moved  out  of  the  interaction 
region  early  in  the  undulator  and  thus  do  not  contribute  to  the  initial  gain  in  the 
undulator.  When  introducing  desynchronism  to  the  SELENE  oscillator  design,  it  is 
imperative  that  the  function  of  the  oscillator  be  maintained.  As  electrons  are  moved 
out  of  the  optical  field,  the  ability  of  the  oscillator  to  bunch  the  electrons  is  reduced. 
Figure  23  shows  the  peak  optical  power  density  at  the  oscillator  mirrors  for  various 
values  of  both  negative  and  positive  desynchronism.  The  optical  power  density  is 
normalized  to  9  MW /cm 2  on  the  vertical  axis,  which  is  equivalent  to  a  peak 
dimensionless  optical  field  of  424.  In  order  to  reduce  the  optical  power  density  to 
50  kW/cm2,  the  amount  of  desynchronism  would  have  to  be  in  excess  of  20%  of  a 
slippage  distance  (0.2  on  the  horizontal  scale),  or  24  optical  wavelengths,  in  advance 
of  the  optical  field.  Typical  values  of  desynchronism  are  on  the  order  of  less  than  1% 


51 


*  SELENE  PEL  Pulse  Evolution  *\ 
;j=i  By^z=5  Oq=Q  c^=O.Oi::'f  <3=0-^f5 

Q=t50  D=0  .  ■  N— 120 


z  8-13  v  13  0  n  200 


Figure  22.  Optical  field  evolution  over  200  passes  for  a  low  current  FEL. 

Parameters:  j  =  1 .0;  az  =  5;  ctg  =  0.01 ;  d  =  0;  Q  =  50;  D  =  0. 

of  a  slippage  distance.  Even  out  to  20%  desynchronism,  the  optical  power  density  is 
in  excess  of  1.2  MW/cm2,  well  above  the  limitations  on  the  mirrors.  Desynchronism, 
by  itself,  will  not  solve  the  problem  of  excessive  gain  in  the  SELENE  oscillator.  Using 
a  desynchronism  of  0.8%  provides  the  lowest  power  output  near  the  typical  (1%) 
desynchronism  level  while  still  maintaining  the  function  of  the  oscillator  in  bunching  the 
electrons. 

Another  way  to  reduce  the  optical  power  in  the  oscillator  is  to  change  the  losses 
in  the  oscillator.  The  major  loss  in  an  oscillator  is  the  semi-transparent  mirror  which 
allows  the  laser  light  to  leave  the  oscillator.  Increasing  the  percentage  of  light 
removed  from  the  oscillator  each  pass  will  result  in  a  lower  steady  state  optical  field  in 
the  oscillator.  For  Q  =  20,  only  5%  of  the  optical  field  power  is  transmitted  through  the 
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Figure  23.  Normalized  power  density  of  the  SELENE  oscillator  as  a  function  of 
desynchronism,  d. 


semi-transparent  mirror.  This  value  of  Q  is  a  typical  value  to  use  for  FEL’s.  Figure  24 
shows  the  result  of  reducing  Q  in  the  SELENE  oscillator.  The  plot  is  normalized  to  a 
peak  power  density  of  19  MW /cm2.  The  optical  power  is  nearly  linearly  proportional  to 
Q  from  Q  -20  down  to  a  value  of  Q  =  2.5.  Figure  25  is  an  expansion  of  Fig.  24 
which  shows  the  relationship  between  optical  power  density  and  Q  for  Q  <  1 .  In  order 
to  reduce  the  optical  power  density  below  50kW/cm2,  Q  must  be  reduced  below 
0.25.  This  means  that  a  transmittance  (1-e_1/Q)  of  98%  is  required  in  the  output 
mirror.  As  can  be  seen  from  Fig.  25,  the  optical  power  drops  rapidly  for  Q  <  0.25. 
This  is  to  be  expected  as  the  gain  in  the  oscillator  is  no  longer  able  to  makeup  for  the 
losses  in  the  cavity  and  the  oscillator  essentially  becomes  a  radiator.  Q  =  5  is  chosen 
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for  a  reasonable  compromise  between  maintaining  low  optical  power  and  keeping 
optical  characteristics  of  the  mirrors  in  a  reasonable  range.  This  results  in  a  20% 
transmittance  for  the  output  mirror. 


Figure  24.  Normalized  power  density  of  the  SELENE  oscillator  as  a  function  of  Q. 

The  original  design  of  SELENE  exploited  the  effect  of  optical  klystrons  in  strong 
optical  fields.  The  dual  dispersive  section  results  in  the  ability  to  change  the  relative 
strength  of  the  dispersive  sections  and  to  allow  for  some  inherent  flexibility  in  the 
design  of  SELENE.  As  previously  shown  in  Fig.  17,  the  single  pass  gain  for  an  optical 
klystron  is  not  much  affected  by  changes  in  dispersive  strength  for  1  <  D  <  4.  This 
can  be  understood  by  looking  at  the  effect  of  the  dispersive  section  on  the  electron 
beam  in  strong  optical  fields.  In  strong  optical  fields,  the  electrons  are  bunched  prior 
to  entering  the  dispersive  section.  If  the  electrons  are  already  bunched  prior  to 
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Figure  25.  Normalized  power  density  of  the  SELENE  oscillator  as  a  function  of  Q  for 
Q<  1. 

entering  the  second  undulator  section,  they  will  be  severely  dispersed  in  the  dispersive 
section.  For  small  values  of  D,  this  dispersive  effect  is  proportional  to  the  dispersive 
strength.  For  larger  values  of  0,  the  dispersion  results  in  a  nearly  random  distribution 
in  the  electrons’  phase;  thus  increasing  the  dispersive  strength  will  result  in  little 
reduction  in  the  single  pass  gain. 

For  SELENE,  the  steady  state  optical  field  strength  is  strong  enough  to  over¬ 
bunch  the  electrons  prior  to  reaching  the  first  dispersive  section.  As  can  be  seen  in 
Fig.  26,  increasing  the  dispersive  strength  of  the  dual  klystrons  will  not  result  in  a 
significant  reduction  in  the  optical  power  at  the  mirrors.  Thus,  with  reasonable 
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selections  of  parameters,  the  excessive  gain  in  the  SELENE  oscillator  cannot  be 
defeated.  The  main  function  of  the  oscillator  is  also  defeated  in  the  current  design. 
Figure  27  shows  the  electron  distribution  in  phase-space  when  they  exit  the  oscillator. 
Instead  of  bunching  the  electrons  for  the  radiator,  the  oscillator  has  essentially 
produced  an  electron  beam  of  the  same  current,  but  with  a  much  larger  energy 
spread.  This  is  a  direct  result  of  the  strong  fields  combined  with  an  optical  klystron  in 
the  oscillator. 


Figure  26.  Normalized  power  density  of  the  SELENE  oscillator  as  a  function  of  the 
dispersive  strength,  D. 


In  order  for  the  oscillator  to  bunch  the  electrons  properly,  the  optical  field  must  be 
reduced,  and  the  dispersive  sections  must  be  tailored  to  the  specific  purpose  of  the 
oscillator.  Using  the  dispersive  section  to  reduce  the  optical  field  may  actually  hinder 
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Figure  27.  Phase-Space  plot  of  the  electron  distribution  exiting  the  SELENE 


oscillator. 


the  ability  of  the  oscillator  to  bunch  the  electrons  for  the  radiator.  The  current  oscillator 
design  does  not  adequately  bunch  the  electrons,  which  leads  to  an  exploration  of 
possible  alternative  designs  for  the  SELENE  oscillator. 
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IV.  ALTERNATE  PROPOSALS  FOR  THE  SELENE  OSCILLATOR 


A.  INTRODUCTION 

In  the  previous  discussion  of  the  oscillator  design  for  SELENE,  some  of  the 
different  effects  of  a  single  dispersive  section  and  single  undulator  design  were 
demonstrated.  Figure  16  shows  that  there  is  a  significant  loss  of  gain  when  only  one 
dispersive  section  (instead  of  two)  is  used  in  weak  optical  fields.  Figure  17 
demonstrates  that  this  may  no  longer  be  the  case  for  strong  field  FEL’s.  In  order  to 
keep  the  flexibility  of  design  inherent  in  the  use  of  an  optical  klystron,  the  two-section 
optical  klystron  will  be  investigated  as  a  first  attempt  to  allow  the  SELENE  oscillator  to 
perform  properly. 


B.  TWO-SECTION  OPTICAL  KLYSTRON  ALTERNATIVE  OSCILLATOR  DESIGN 

Figure  21  shov/s  that  just  using  a  single  dispersive  section  in  the  same  undulator 
design  will  still  produce  excessive  gain  and  strong  optical  fields  (a  »  840).  Thus,  the 
need  to  reduce  gain  leads  to  the  necessity  for  a  shorter  undulator.  Reviewing  the 
previous  discussions  of  the  relevant  parameters  leads  to  the  proportionalities. 


oG  <*  N  , 
oQ<*  N  , 


oc 


N  ’ 


j  «  A/3  , 

laU«  J-M.lSl.Js  . 

zo 


(4.1) 


In  order  to  determine  the  length  of  undulator  desired  for  the  two-section  oscillator,  the 
best  parameter  to  use  is  the  dimensionless  current  density.  With  low  current 


59 


conditions,  the  gain  in  a  single  pass  undulator  is  given  by  G  =  jD/4.  Thus  with  j  =  1 , 
the  gain  in  a  single  pass  is  D/4.  With  proper  choices  for  desynchronism  and 
resonator  losses  small  values  of  D  may  enable  the  oscillator  to  remain  in  the  weak 
field  regime.  Using  Eq.  4.1,  the  parameters  used  to  describe  the  electron  beam 
distribution  and  slippage  distance  for  j  =  1  are  N  =  26,  oG  =  0.15,  a9  =  0.07  ,  and 
oz  =  17.5.  Again  assuming  the  mirrors  are  designed  to  produce  a  Rayleigh  length 
equal  to  the  length  of  one  undulator  section,  z0  «  1.1  m.  This  will  result  in  a  maximum 
dimensionless  optical  field  strength  of  la lmax  =  4.6  in  order  not  to  exceed  50  kW/cm2 
at  the  mirrors,  keeping  the  mirror  separation  distance  at  Sm  =  79  m. 

Figure  28  is  a  plot  of  the  normalized  peak  power  density  after  1000  electron 
pulses  as  a  function  of  desynchronism.  The  peak  power  density  is  normalized  to 
211  kW/cm2.  As  before,  negative  desynchronism  results  in  the  greatest  decrease  in 
optical  power  for  small  changes  in  desynchronism  (steepest  slope).  However, 
negative  desynchronism  will  prevent  a  larger  number  of  electrons  from  interacting  with 
the  optical  field.  For  d  >  0,  the  peak  power  density  decreases  as  desynchronism  is 
increased.  Thus  a  small  value  of  positive  desynchronism  of  d  =  0.02  is  best,  in  this 
case,  to  ensure  the  largest  number  of  electrons  continue  to  interact  with  the  optical 
field  while  allowing  for  a  small  reduction  in  the  steady  state  optical  field  strength. 

After  determining  the  desynchronism  to  be  used  in  this  two-section  oscillator,  the 
next  step  is  to  determine  the  dispersive  strength  which  will  optimize  the  performance 
of  the  oscillator  and  reduce  the  optical  field.  Figure  29  is  a  dual  plot  showing  the 
normalized  peak  power  density  and  the  normalized  peak  dimensionless  optical  field 
strength  as  a  function  of  dispersive  strength,  D.  The  peak  power  density  is 
normalized  to  284  kW/cm2  and  the  peak  dimensionless  optical  field  strength  is 
normalized  to  I  a  I  =  1 1 .  The  peak  power  is  seen  to  be  reduced  as  the  dispersive 
strength  is  increased.  The  peak  dimensionless  optical  field  strength  follows  a  similar 
reduction  until  D  >  0.3.  At  this  dispersive  strength,  the  optical  field  tends  to  be 
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Figure  28.  Normalized  peak  power  density  as  a  function  of  desynchronism  for  the 
two-section  optical  klystron  configured  for  y  =  1 . 

reduced  much  more  slowly  than  the  peak  optical  power.  This  can  be  explained  using 
Figs.  30  and  31 . 

Figure  30  shows  the  optical  pulse  evolution  over  3000  passes.  It  can  clearly  be 
seen  from  the  upper  left-hand  plot,  that  the  optical  pulse  shape  is  smooth  and  nearly 
symmetrical.  An  increase  in  the  peak  optical  field  strength  with  this  optical  envelope 
shape  will  be  accompanied  by  an  increase  in  the  average  power  in  the  oscillator. 
However,  Fig.  31  shows  the  same  oscillator  with  a  dispersive  strength  of  D  =  0.4.  In 
this  figure,  the  reason  for  the  discrepancy  in  Fig.  29  becomes  apparent.  The  stronger 
dispersive  strength  in  the  presence  of  strong  optical  fields  tends  to  produce  many 
small  areas  of  bunching  along  the  optical  pulse  and  therefore,  the  optical  field  is 
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Figure  29.  Normalized  peak  power  density  and  optical  field  strength  as  a  function  of 
dispersive  strength  for  the  two-section  optical  klystron  with  j  =  1 . 

amplified  in  a  non-uniform  manner.  This  causes  the  spikes  in  the  optical  field 
envelope.  This  can  be  avoided  by  using  lower  values  of  dispersive  strength,  or 
lowering  the  optical  field.  For  this  situation,  the  aberrations  in  the  optical  field 
envelope  lead  to  the  decision  to  choose  D  =  0.3  for  a  dispersive  strength.  This  will 
produce  a  decrease  in  power  density  to  ~  65%  of  the  steady  state  power  density 
attained  without  the  dispersive  section,  or  to  185  kW/cm2.  Thus,  a  decrease  in  the 
dimensionless  optical  field  strength  by  a  factor  of  two  is  still  necessary  in  order  to 
prevent  damage  to  the  mirrors  in  the  oscillator. 

Figure  32  shows  the  normalized  peak  power  density  as  a  function  of  Q.  The 
power  density  is  normalized  to  177  kW/cm2.  The  optical  field  intensity  is  reduced 
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Figure  30.  Optical  pulse  evolution  during  initial  3000  electron  pulses.  D  -  0.3  for 
the  two-section  optical  klystron. 

rapidly  as  the  resonator  losses  are  increased.  For  Q  <  6.6,  the  undulator  gain  is 
insufficient  to  make  up  for  the  resonator  losses.  If  the  FEL  oscillator  output  mirror  is 
chosen  such  that  Q  ~  7,  then  the  dimensionless  optical  field  strength  is  low  enough  to 
prevent  damage  to  the  mirrors  in  the  oscillator.  This  seems  to  be  a  reasonable 
proposal  to  allow  the  oscillator  to  work  but  not  damage  the  mirrors  in  the  oscillator. 
However,  the  function  of  the  oscillator  must  be  considered. 

Figure  33  shows  the  phase-space  plot  for  the  electrons  at  the  end  of  the 
oscillator.  The  electrons  plotted  in  light  gray  are  the  positions  of  the  electrons  just 
prior  to  the  dispersive  section.  The  electrons  plotted  in  dark  gray  are  the  positions  of 
the  electrons  as  they  exit  the  undulator.  As  can  be  seen  in  comparison  with  Fig.  1 1 , 
the  dispersive  section  acts  to  bunch  the  electrons.  Since  the  initial  optical  field  is  not 
strong  enough  to  have  over-bunched  the  electrons  prior  to  the  dispersive  section,  the 
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Figure  31.  Optical  pulse  evolution  during  the  initial  3000  electron  pulses  with 
D  =  0.4  for  the  two-section  optical  klystron. 

dispersive  section  in  the  oscillator  bunches  the  electrons  and  the  second  undulator 
section  continues  to  interact  with  the  optical  field  and  the  results  are  seen  in  Fig  33. 

At  first  glance,  the  objective  of  the  oscillator  seems  to  have  been  met.  The 
electrons  are  bunched  at  the  end  of  the  oscillator  and  the  mirror  damage  threshold  is 
not  approached.  However,  there  is  still  the  small  dispersive  section  between  the 
oscillator  and  the  radiator  to  consider.  As  previously  discussed  and  shown  in  Fig  13, 
having  the  electrons  bunched  prior  to  entering  a  dispersive  section  results  in  the 
electrons  being  dispersed  instead  of  bunched.  This  will  result  in  reduced  efficiency  in 
the  radiator  in  extracting  energy  from  the  electron  beam.  If,  however,  upon  exiting  the 
oscillator,  the  electron  beam  could  look  similar  to  the  electrons  plotted  in  light  gray  in 
Fig  33,  then  the  dispersive  section  between  the  oscillator  and  the  radiator  will  bunch 
the  electrons  and  the  radiator  efficiency  will  be  increased  dramatically.  This  may  be 
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accomplished  by  removing  the  dispersive  section  in  the  oscillator  and  reducing  the 
dimensionless  optical  field  strength  to  the  weak  field  regime. 


C.  SINGLE  UNDULATOR  ALTERNATIVE  OSCILLATOR  DESIGN 

As  stated  before,  removing  the  dispersive  section  from  the  previous  oscillator 
design  is  not  sufficient  to  provide  the  electron  beam  conditioning  desired  for  the 
radiator.  The  presense  of  strong  optical  fields  will  result  in  the  electrons  being 
bunched  prior  to  leaving  the  oscillator.  The  condition  desired  is  similar  to  Fig.  4.  The 
sinusoidal  shape  of  the  phase-space  plot  is  the  most  desirable  condition  for  the 
electron  beam  to  exit  the  oscillator.  This  particular  proposed  design  will  result  in  the 
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Figure  33.  Phase-Space  plot  for  the  two-section  optical  klystron. 


electron  beam  being  bunched  prior  to  entering  the  radiator  instead  of  at  the  exit  of  the 
oscillator.  In  order  to  reduce  the  optical  field  to  the  weak  field  regime,  the  previous 
design,  without  the  dispersive  section  and  a  lower  Q  could  be  used.  However,  using 
a  lower  value  of  Q  will  tend  to  make  the  oscillator  unstable.  As  can  be  seen  in  Fig. 
32,  the  output  power  is  sensitive  to  changes  in  resonator  losses  near  the  value  which 
will  produce  the  required  weak  optical  field.  In  order  to  save  the  cost  of  constructing  a 
larger  undulator  for  the  oscillator,  and  to  try  and  find  a  less  sensitive  solution,  the 
undulator  may  be  shortened. 

An  undulator  length  of  2  m  is  a  reasonable  solution  to  this  problem.  With  this 
design,  and  using  Eq.  4.1 ,  the  electron  beam  parameters  are  N  =  22,  j  =  0.6, 
aG  =  0.12,  a9  =  0.06,  oz  -21,  and  la  lmax  =  2.0.  Again,  desynchronism  and  resonator 
losses  will  play  an  important  part  in  this  design.  Figures  34  and  35  are  indications  of 
why  these  two  factors  are  so  important.  In  Fig.  34,  the  combination  of  resonator 
losses  and  slippage  results  in  the  FEL  failing  to  sustain  itself.  The  gain  is  small 
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compared  to  the  original  design  of  the  SELENE  oscillator.  The  interaction  distance  is 
also  short,  thus  making  the  slippage  of  the  electron  beam  along  the  optical  beam  an 
important  feature.  In  the  beginning  of  a  simulation  of  this  type,  the  optical  field  is 
seeded  with  a  weak  field  of  strength  a0  =  0.001 .  Early  in  the  evolution,  the  electron 
micropulse  interacts  with  the  optical  field  resulting  in  amplification  of  the  field.  Due  to 
the  slippage  of  the  electron  pulse,  the  optical  field  tends  to  be  amplified  along  the 
trailing  edge.  Due  to  resonator  losses,  the  overall  gain  for  the  oscillator  is  small.  As 
seen  in  the  middle-left  plot,  this  makes  the  optical  field  "walk"  backwards.  After 
*  1 000  electron  pulses,  the  optical  field  and  the  electron  beam  are  no  longer 
interacting  sufficiently  to  produce  gain  in  the  optical  field. 


Figure  34.  Optical  pulse  evolution  over  5000  electron  pulses  for  the  single  undulator 
segment  design  with  Q  =  20  and  d  =  0. 

In  Fig.  35,  the  resonator  losses  are  reduced  so  that  Q  =  25.  This  is  sufficient  to 
enable  the  electron  beam  to  interact  and  produce  a  stable  optical  field  over  many 


67 


passes.  However,  the  shape  of  the  optical  field  is  not  optimum.  This  "shark-fin" 
shape  will  produce  non-uniform  bunching  in  the  electron  beam.  Thus,  the  need  for 
desynchronism.  Due  to  the  nature  of  the  shape  of  the  optical  field,  as  desynchronism 
is  increased,  the  magnitude  of  the  peak  optical  field  will  have  a  small  initial  increase 
since  a  larger  fraction  of  the  electrons  will  spend  time  interacting  with  the  peak  optical 
field.  As  this  happens,  the  optical  wave  shape  will  expand.  Because  of  these  two 
phenomena,  as  desynchronism  is  increased  near  d  =  0,  the  peak  optical  field  will 
remain  relatively  constant  while  the  total  power  in  the  optical  field  increases 
significantly. 
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Figure  35.  Optical  pulse  evolution  over  8000  electron  pulses  for  the  single  undulator 
segment  design  with  Q  =  25  and  d  =  0. 


This  phenomenon  will  continue  until  a  symmetric  waveform  is  approached  as  in 
Fig.  36.  At  d  =  0.02,  the  shape  of  the  optical  pulse  is  symmetric  and  uniform. 
Increasing  desynchronism  beyond  this  point  will  provide  the  same  results  as  discussed 
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in  previous  sections.  As  desynchronism  is  increased,  the  peak  optical  field  will 
decrease  and  the  total  power  in  the  optical  field  will  also  decrease.  Note  that  negative 
desynchronism  in  this  case  would  be  counter-productive  as  it  would  remove  the 
electron  pulse  from  the  optical  field  envelope  much  sooner  than  without 
desynchronism.  This  will  not  produce  a  stable  optical  field  and  is  not  a  choice  for  this 
system.  Starting  with  the  uniform  optical  field  envelope  from  Fig.  36,  Fig.  37  can  be 
constructed.  Figure  37  is  a  plot  of  peak  optical  field  strength  as  a  function  of  Q.  With 
d  =  0.02,  the  shape  of  the  optical  field  pulse  remains  consistent  as  Q  is  decreased. 
At  Q  =  13.8  the  peak  optical  field  strength  is  reduced  below  la  I  =  2  which  will  prevent 
damage  to  the  conservative  mirror  elements.  This  also  is  in  the  weak  field  regime. 


Figure  36.  Optical  pulse  evolution  over  8000  electron  pulses  for  the  single  undulator 
segment  design  with  Q  =  25  and  d  =  0.02. 

Figure  38  is  the  phase-space  plot  of  the  electrons  exiting  the  oscillator  cavity  with 
an  optical  field  strength  of  a  =  2.0.  This  is  the  type  of  sinusoidal  shape  desired  in 
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Q 

Figure  37.  Dimensionless  optical  field  strength  as  a  function  of  Q  for  the  single 
undulator  segment  design. 

order  to  take  advantage  of  the  dispersive  section  between  the  oscillator  and  the 
radiator.  However,  it  will  be  important  to  choose  the  dispersive  section  strength 
properly  in  order  to  continue  the  bunching  process  without  overbunching.  With  the  use 
of  an  achromatic  bending  system,  the  electron  beam  can  be  introduced  into  the 
radiator  without  any  further  bunching.  A  compromise  between  the  achromatic  bending 
and  the  "current  designed  dispersive  strength  is  required  for  optimum  performance  of 
the  radiator. 

For  further  consideration,  the  flexibility  of  this  system  can  be  expanded  upon. 
With  increasing  advances  in  the  design  of  optical  devices,  mirrors  which  can  withstand 
greater  than  200kW/cm2  of  incident  optical  power  density  may  be  available  by  the 
time  the  SELENE  project  is  finalized.  This  will  expand  the  region  of  acceptable 
resonator  loss  to  13.6  <  Q  £  15.  With  the  use  of  a  feedback  and  control  system,  the 
desynchronism  of  the  oscillator  may  be  dynamically  controlled.  With  the  main  purpose 
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Figure  38.  Phase-Space  plot  for  the  electrons  exiting  the  single  undulator  segment 
oscillator. 

of  the  oscillator  being  to  maintain  weak  field  electron  beam  bunching,  a  feedback  loop 
may  be  implemented  to  dynamically  control,  the  desynchronism  to  maintain  the  optical 
field  in  the  weak  field  regime. 

This  feedback  system  can  also  be  used  in  conjunction  with  a  variable  dispersive 
section  system  between  the  oscillator  and  the  radiator.  This  design  originally 
proposed  for  the  dispersive  sections  inside  the  oscillator  can  be  used  to  assist  in 
dynamically  controlling  electron  bunching.  As  the  optical  field  strength  decreases,  the 
strength  of  the  dispersive  section  can  be  increased  to  bunch  the  electrons  prior  to 
entering  the  radiator.  As  the  optical  field  increases,  the  dispersive  section  can  be 
made  achromatic,  in  other  words,  non-dispersive.  As  the  dimensionless  optical  field 
strength  reaches  a  =  4,  the  dispersive  section  should  be  fully  achromatic.  With  full 
use  of  an  integrated  feedback  system,  the  output  from  the  radiator  can  be  optimized 
dynamically  and  thus  the  available  power  to  the  satellites  will  be  maximized. 
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V.  CONCLUSIONS 


With  the  current  explosion  in  electronics  technology,  the  United  States  military 
services  are  developing  a  large  arsenal  of  data  gathering  sensors  and  equipment. 
This  data,  in  most  cases  is  gathered  in  the  field  and  needs  to  be  transmitted  to  a 
coordinating  unit  for  processing.  Advances  in  optical  devices  and  sensors  are 
increasing  the  range  and  reliability  of  optical  sensors  in  various  ranges  of  frequencies 
from  microwave  to  visual  imaging.  The  number  of  sensors  in  the  field  which  use  video 
technology  for  relaying  information  is  increasing  rapidly.  The  communications 
requirement  for  the  military  services  is  expanding.  The  days  of  the  walkie-talkie 
communications  are  numbered.  All  of  these  communications  needs  must  be  met.  The 
bandwidth  required  to  transmit  video  images  on  a  real-time  basis  is  an  order  of 
magnitude  larger  than  that  required  to  transmit  voice.  Navigation  and  precise 
positioning  information  is  also  necessary  for  current  targeting  data  to  prevent  collateral 
damage.  With  the  rapid  proliferation  of  cellular  phone  systems,  the  commercial 
communications  world  is  also  expanding  it’s  needs  to  provide  world-wide  continuous 
coverage  for  customers.  All  of  these  concerns  have  at  least  one  thing  in  common. 
They  all  benefit  from  the  use  of  satellites. 

In  the  near  future,  a  global  telecommunications  system  will  be  installed  using  a 
series  of  satellites.  The  number  of  people  requiring  continuous  Global  Positioning 
SystenT(GPS)  service  will  rapidly  expand  as  different  applications  are  introduced  into 
society.  Currently,  even  some  novice  of  fishermen  have  a  GPS  system  on  board 
mounted  next  to  the  fish-finder.  In  the  near  future,  it  can  be  imagined  that  cars  can  be 
programmed  with  a  travel  route  and  with  the  use  of  GPS  and  other  appropriate 
sensors,  a  long  arduous  road  trip  may  be  made  without  error,  and  possibly  without 
manual  intervention.  The  newest  television  systems  are  using  satellite  technology  to 
provide  interactive  programming  directly  to  a  small  communications  satellite  dish. 
Because  of  the  increasing  reliance  on  satellite  systems,  and  the  increasing  functional 
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requirement  of  these  satellites,  any  method  which  will  lower  the  orbital  placement  cost 
of  new  satellites,  or  extend  the  lifetimes  of  current  satellites  will  be  welcome  by  both 
industrial  interests  and  by  consumers.  The  military  services,  which  use  both  special 
military  satellites  and  commercial  units,  will  benefit  greatly  from  such  an  improvement. 
This  advanced  technology  is  the  purpose  of  SELENE. 

The  United  States  and  Russian  governments  have  agreed  to  conduct  a  joint 
development  of  a  ground-based,  high-average-power  laser  system  to  provide  electrical 
power  to  satellites  in  orbit.  The  design  will  use  a  Free  Electron  Laser  (FEL)  system. 
This  FEL  system  consists  of  two  main  components  in  addition  to  the  electron  beam 
source:  an  oscillator,  and  a  radiator.  The  electron  beam  will  be  provided  by  a  51  MeV 
race-track  microtron  recuperator.  The  electron  beam  will  be  injected  into  the  oscillator 
which  will  bunch  the  electrons.  The  electrons  will  then  enter  the  radiator  where  they 
will  produce  a  coherent  optical  field  by  spontaneous  emission.  After  the  field  is 
developed,  it  is  further  amplified  by  the  interaction  in  the  undulator  between  the 
electrons  and  the  optical  field.  The  resulting  laser  light  is  then  transmitted  through  a 
two-miie-long  vacuum  chamber  to  allow  it  to  expand  before  reaching  the  beam 
director.  The  beam  director  will  reflect  the  laser  beam  to  the  satellites  using  adaptive 
optics  techniques  to  compensate  for  atmospheric  disturbances. 

In  order  to  provide  high  power  from  the  radiator,  the  electron  beam  current  must 
be  large.  In  the  oscillator,  this  causes  several  problems.  Mirror  damage  is  a  concern 
since  strong  optical  fields  could  be  developed  in  the  oscillator.  These  strong  optical 
fields  tend  to  destroy  the  bunching  of  the  electrons  as  they  travel  down  the  oscillator. 
The  original  three-section  optical  klystron  was  used  to  try  and  reduce  the  strong  optical 
fields  in  the  oscillator.  Chapter  III  discusses  the  proposed  system  for  the  SELENE 
oscillator  and  the  limitations  of  the  system.  With  current  mirror  technology,  it  is  difficult 
to  develop  a  stable,  optical  field  using  that  system  that  will  not  exceed  the  damage 
threshold  of  the  mirrors.  Even  if  the  mirrors  could  be  designed  to  withstand  the 
intense  optical  field  produced  by  this  oscillator,  the  purpose  of  the  oscillator  would  not 
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be  served  for  the  following  reason.  Figure  13  shows  the  result  of  using  dispersive 
sections  in  strong  optical  fields.  The  electrons  are  bunched  prior  to  reaching  the 
dispersive  section.  When  the  electrons  exit  the  dispersive  section,  the  electrons  have 
been  scattered  instead  of  being  bunched.  SELENE  has  not  just  one,  but  two, 
dispersive  sections  operating  in  strong  fields.  The  result  is  shown  in  Fig.  27.  The 
electrons  leaving  the  oscillator  are  not  bunched.  Instead,  this  combination  of 
dispersion  and  strong  optical  fields  produces  a  relatively  uniform  electron  distribution  in 
C  with  a  much  larger  energy  spread  than  the  original  electron  beam.  This  will  inhibit 
the  creation  of  coherent  radiation  in  the  radiator. 

As  an  alternative  to  the  SELENE  proposal,  a  two-section  optical  klystron  was 
studied.  The  two  undulator  sections  were  shortened  in  order  to  reduce  the 
dimensionless  current  density  to  y  =  1.  The  resulting  investigation  showed  that  by 
proper  selection  of  dispersive  strengths,  desynchronism,  and  resonator  losses,  the 
goal  of  weak  field  bunching  in  the  oscillator  could  be  met.  However,  this  solution  is 
very  sensitive  to  resonator  losses.  The  other  concern  is  that  even  under  the  weak 
field  conditions,  the  electrons  are  bunched  when  they  exit  the  oscillator.  The 
dispersive  section  between  the  oscillator  and  the  radiator  will  again  disperse  the 
electrons  instead  of  bunching  them. 

Another  alternative  to  the  SELENE  proposal  is  to  remove  the  dispersive  section 
and  shorten  the  undulator.  At  first,  this  configuration  seems  to  require  a  large  Q  in 
order  to  produce  a  stable  optical  field.  Once  a  value  for  resonator  losses  was  found 
which  allowed  the  production  of  an  optical  field,  the  shape  of  the  optical  pulse  became 
an  issue.  Selecting  a  desynchronism  value  which  optimized  the  shape  of  the  optical 
field  resulted  in  d  =  0.02.  With  this  desynchronism,  the  resonator  losses  were 
increased  until  a  stable,  weak  optical  field  was  formed.  As  Q  was  reduced,  the 
number  of  passes  required  to  form  the  optical  field  grew,  but  the  optical  field  remained 
stable  as  the  steady  state  dimensionless  optical  field  strength  was  reduced  to  a  =  1 .3, 
well  within  the  weak  field  range.  For  13.5  £  Q  £14.3,  the  single  section  oscillator  is 
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able  to  produce  a  stable  weak  optical  field.  In  Fig.  39,  the  optical  pulse  evolution  over 
50,000  electron  pulses  is  shown  for  Q  =  14.  The  optical  field  stability  in  the  weak  field 
range  allows  for  dynamic  control  of  desynchronism  and  the  strength  of  the  dispersive 
section  between  the  oscillator  and  the  radiator. 
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Figure  39.  Optical  pulse  evolution  during  50000  electron  pulses  for  the  one-piece 
undulator. 


Figure  40  demonstrates  the  usefulness  of  the  dispersive  section  between  the 
oscillator  and  the  radiator.  The  left-hand  side  is  the  phase-space  plot  of  the  electrons 
as  they  exit  the  oscillator.  The  right-hand  side  is  the  phase-space  plot  of  the  electrons 
as  they  enter  the  radiator  after  passing  through  a  dispersive  section  with  D  =  0.25. 
Note  that  the  effect  of  the  dispersive  section  in  weak  optical  fields  is  to  bunch  the 
electrons.  Because  the  radiator  does  not  have  an  initial  optical  field,  the  relative 
phase  at  which  the  electrons  are  bunched  is  not  important,  just  that  they  are  bunched 
in  phase  as  they  begin  to  produce  coherent  radiation  by  spontaneous  emissions. 
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a)  Exit  of  the  oscillator.  b)  Entrance  to  the  radiator 


Figure  40.  Phase-Space  plot  of  the  electrons  as  they  a)  exit  the  oscillator  and  b) 
enter  the  radiator  of  SELENE.  The  strength  of  the  dispersive  section  is 
D  =  0.25. 

Using  a  variable  strength  for  this  dispersive  section,  from  an  achromatic  bend  for 
a  >  4  to  D  »  1  for  weak  fields  of  a  ~  1.5  will  allow  for  dynamic  control  to  maximize 
bunching  of  the  electrons  prior  to  entering  the  radiator.  This  design  meets  the 
requirements  of  bunching  the  electrons  for  the  radiator  and  reducing  the  optical  field 
strength  below  the  damage  threshold  of  current  mirror  technology. 
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